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SECTION  I 


INTRODUCTION 

This  investigation  is  concerned  with  the  prediction  of  the  nonlinear 
response  of  reinforced  concrete  structures,  including  member  failures  and 
structural  collapse,  to  static  and  dynamic  loads.  The  computer  program 
SINGER,  the  product  of  this  investigation,  provides  the  tool  for  this 
prediction.  This  report  describes  the  mathematical  models  and  the  solution 
process  which  form  the  basis  of  SINGER. 

1.1  BACKGROUND 

Since  it  was  desired  to  represent  the  structure  by  a  discrete  model 
composed  of  "gross  elements"  it  was  natural  to  select  the  finite  element 
method  to  model  the  structure.  However,  the  selection  of  the  solution 
process  represented  a  pivotal  decision.  Two  methods  were  given  serious 
consideration:  the  step-by-step  (STEP)  approach,  an  equilibrium  approach 

in  which  the  structure  is  represented  by  a  stiffness  matrix;  and  the 
minimization  (MIN)  approach,  an  energy  approach  in  which  the  structure 
is  characterized  by  a  work  function.  In  both  approaches,  the  solution 
process  initiates  at  a  point  where  the  state  of  the  system  is  known 
and  proceeds  along  discrete  points  of  the  equilibrium  path  (motion)  of 
the  system. 

The  STEP  approach  has  been  used  extensively  in  the  analysis  of  non¬ 
linear  structures  and  is  well  documented  [e.g.,  14,17]*.  The  central  idea 
of  the  STEP  approach  is  contained  in  Newton's  method  of  successive  ap¬ 
proximations  to  a  real  root  [15],  It  is  illustrated  in  Figure  la, 
which  depicts  the  nonlinear  equilibrium  path  of  a  one-degree-of-freedom 

* 

Numbers  in  brackets  designate  references 
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during  the  entire  iterative  process  (see  Figure  lb)  or  by  combining  the 
constant  and  variable  stiffness  coefficients  in  the  solution  process  [17]. 
The  extension  of  Newton's  method  to  a  system  with  multi-degrees-of- 

freedom  is  known  as  the  Newton-Raphson  method.  On  the  basis  of  the  finite 
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element  method,  the  governing  equations  of  equilibrium  can  be  expressed 
in  the  form  [17] 

p  *  /B^odV  (1.8) 

V 

where 

c  =  Bx  (1.9) 

and 

6e  *  B6x  (1.10) 

In  Equations  1.8,  1.9,  and  1.10,  p  and  x  represent  the  external  generalized 
force  and  displacement  vectors,  respectively;  c  and  e  denote  the  stress  and 
strain  vectors,  respectively  (the  constitutive  laws  may  be  nonlinear); 

B  is  th'  compatibility  matrix  which  may  depend  on  x,  in  which  case 
B  +  B;  6  signifies  a  virtual  variation;  and  V  denotes  the  volume  of 
the  system.  Again  the  condition  of  equilibrium  for  a  specific  force 
vector  p  is 

Ap  =  p  -  /B^adV  *  0  (1.11) 

V 

where  Ap  is  the  unbalanced  force  vector.  Analogous  to  the  Newton  process, 

the  correction  to  the  n^  trial  solution,  x  ,  is 

n 

Ax  =  K_-1Ap  (1.12) 

n  l  n 

s  t 

pnd  the  n  +  1  trial  solution  is  defined  by 

x  , ,  =  x  +  Ax  (1.13) 

n+1  n  n 

The  tangent  stiffness  matrix  K^,  in  Equation  1.12  is  obtained  by  forming 
a  virtual  variation  of  Equation  1.8  with  respect  to  x;  the  result  can  be 
expressed  in  the  form 

<5p  =  I^ox  •  (1.14) 
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The  unbalanced  f 01 ce  vector  c o  responding  to  any  trial  solution  is 
evaluated  on  the  basis  of  Equation  1.11.  The  solution  process  is 
continued  until  r.ha  unbalanced  forces  are  sufficiently  small.  The 
modifications  of  the  Nc-wton  process  are  also  employed  in  the  Newton- 
Raphson  process. 

The  MIN'  anproach  is  based  on  the  oroperty  that  the  work  function 
[7]  of  the  systen  assumes  a  relative  minimum  at  a  stable  equilibrium 
state.  Accordingly,  a  desired  equilibrium  state  is  found  by 
minimization  of  the.  work  function.  Function  minimization  is  accom¬ 
plished  via  nonlinear  programming  techniques.  The  MIN  approach 
has  been  employed  successfully  in  the  analysis  of  nonlinear  structures 
[e.g. ,  2,  5,  9] . 

The  Mill  process,  which  is  discussed  in  more  detail  in  section  3, 

is  illustrated  for  a  two-degree-of -freedom  system  in  Figure  2.  The 

work  function  W  is  represented  by  level  curves.  Function  minimization 

is  based  on  a  modification  of  Davidon's  method  [13].  The  search  for 

the  desired  equilibrium  state  x  corresponding  to  the  applied  load 

vector  p  Initiates  at  x  in  the  direction  d. .  The  first  trial  solution 

o  i 

x  is  obtained  by  minimizing  the  function  W  along  the  direction  d^.  A 
new  search  direction  d9  is  established,  and  the  relative  minimum  of 
W  with  respect  to  d,,  is  found  to  be  x2  (the  search  directions  are  defined 
by  transformations  of  the  gradients  of  the  work  function  [13]).  The 
iterative  process  is  continued  until  the  components  of  the  gradient 
of  the  work  function,  which  correspond  to  the  unbalanced  forces,  are 
sufficiently  small. 

It  was  decided  that  both  the  STEP  and  MIN  approach  provide  a 
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satisfactory  basis  for  the  proposed  modeling  and  solution  process.  The 
final  selection  of  the  MIN  approach  was  strongly  influenced  by  the  fol¬ 
lowing  factors: 

Ir.  the  MIN  process,  the  search  for  an  equilibrium  state  is  always 
based  on  the  actual  state  of  the  system  corresponding  to  an  assumed 
displacement  configuration.  The  STEP  approach  is  a  quasi-linear 
approach  in  which  every  trial  solution  is  based  on  the  stiffness 
properties  of  the  system  at  the  beginning  of  the  iteration  (Figure  1) . 
Hence,  in  the  MIN  approach,  decisions  are  always  based  on  the  actual 
state  of  the  system. 

Depending  on  the  choice  of  the  minimization  algorithm,  substantial 
storage  space  savings  can  be  achieved  with  the  MIN  process  since 
the  structure  is  represented  by  a  scalar  function.  However,  the 
Davidon  algorithm  [13]  does  require  storage  space  comparable  to 
the  STEP  process. 

The  computation  and  "inversion"  of  the  tangent  stiffness  matrix  in 
the  Newton-Raphson  process  (Equation  1.12)  requires  a  significant 
amount  of  computational  effort  at  each  cycle.  The  alternate  approach 
of  a  constant  stiffness  matrix  converges  only  for  certain  types  of 
nonlinearitres  [14].  Hence,  at  least  a  combination  of  the  constant 
and  variable  stiffness  matrices  is  required. 

1 . 2  PURPOSE  AND  SCOPE 

The  function  of  the  computer  program  SINGER  is  to  predict  the  be¬ 
havior  of  plane  skeletal  reinforced  concrete  structures  in  their  environ- 
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ments.  Of  particular  interest  is  the  nonlinear  transient  response  including 
the  possibility  of  element  failures  and  structural  collapse. 

SINGER  is  intended  to  serve  as  a  tool  fcr  the  improvement  and  devel¬ 
opment  of  techniques  for  the  assessment  of  existing  protective  structures, 
the  design  of  new  systems,  and  the  development  of  motion  environment 
criteria  for  internal  systems  of  protective  structures. 

1.3  METHODOLOGY 

The  prediction  of  the  performance  of  the  structure  in  its  environ¬ 
ment  is  based  on  the  response  of  a  mathematical  model  of  the  structure  to 
actions,  which  simulate  the  environment.  The  analysis  process  comprises 
three  principal  tasks: 

1.  The  formulations  of  actions,  the  mathematical  models  of  the 
environment. 

2.  The  development  of  a  mathematical  model  of  the  structure. 

3.  The  formulation  of  the  solution  process. 

The  actions  consist  of  the  self-weight  of  the  structure,  distributed 
and  concentrated  static  and  dynamic  loads,  inertia  forces,  and  support 
motions. 

The  structure  is  represented  by  an  assemblage  of  discrete  line  ele¬ 
ments  and  springs  interconnected  at  a  finite  number  of  points.  The  line 
elements  are  models  of  straight,  prismatic,  reinforced  concrete  members 
whose  longitudinal  plane  of  symmetry  corresponds  to  the  plane  of  loading. 

The  line  element  is  discretized  via  the  finite  element  method;  the  inter¬ 
nal  energy,  which  characterizes  the  state  of  the  element,  is  a  function 
of  the  element  distortion  components  (three  relative  end-displacements 
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and  one  relative  internal- displacement) .  Springs  represent  models  of 
joints  with  partial  releases.  A  concentrated  mass  is  assigned  to  each 
degree  of  freedom  of  the  assemblage.  Energy  dissipation  resulting  from 
inelastic  behavior  accounts  for  structural  damping. 

The  line  elements  admit  geometric  and  physical  nonlinearities. 
Geometric  nonlinear i ties  are  induced  by  the  coupling  of  flexural  and 
axial  distortions  and  the  formulation  of  equilibrium  for  the  deformed 
state  of  the  assemblage.  Physical  nonlinearities  are  caused  by  non¬ 
linear  constitutive  (stress-strain)  laws.  The  springs  are  assumed  to 
behave  linearly. 

The.  behavior  of  the  element  is  modeled  up  to  the  limit  of  contin¬ 
uous  change  of  state,  defined  as  fracture  (e.g., crushing  of  the  compres¬ 
sion  block  constitutes  element  failure;  however,  minor  discontinuties 
such  as  spalling  of  the  concrete  cover  are  modeled). 

In  the  linear  domain,  the  state  of  the  system  is  completely  defined 
by  the  generalized  coordinates  which  consist  of  nodal  displacements, 
relative  internal  element-displacements,  and  relative  release-displace¬ 
ments.  In  the  nonlinear  range,  the  generalized  coordinates  must  be 
related  to  the  motion  (equilibrium  path)  of  the  system  to  define  the 
state  of  the  system.  The  origin  of  the  generalized  coordinates  cor¬ 
responds  to  the  unstrained  state  of  the  system,  termed  the  initial  state. 

The  response  of  the  system  to  dynamic  actions  is  determined  at  a 
discrete  number  of  points  in  time.  The  solution  process  is  a  closed 
iterative  process  within  two  successive  points  in  time,  the  time  step. 

The  time  function  of  each  generalized  coordinate  is  approximated 
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over  the  time  step  by  a  finite  power  series  whose  coefficients  are  ex¬ 
pressed  in  terns  of  three  known  initial  conditions,  the  displacement, 
velocity,  and  acceleration  at  the  beginning  of  the  time  step,  and  one 
unknown  end  condition,  the  displacement  at  the  end  of  the  time  step. 

This  representation  of  the  time  function  permits  one  to  express  the  in¬ 
ertia  forces  at  the  end  of  the  time  step  in  terms  of  the  unknown  dis¬ 
placements.  Consequently,  the  state  of  the  system  at  the  end  of  the  time 
step  can  be  completely  defined  in  terms  of  the  corresponding  generalized 
coordinates.  For  this  purpose  a  work  function  is  introduced,  a  scalar 
function  of  the  generalized  coordinates,  which  contains  implicitly  all 
the  forces  acting  on  the  system  (applied,  inertia,  internal).  The  de¬ 
sired  system  configuration  at  the  end  of  the  time  step  is  obtained  by 
minimization  of  the  work  function,  which  assumes  a  relative  minimum  at 
the  dynamic  equilibrium  state.  The  minimization  process  is  a  search 
process  in  which  a  system  configoration  is  assumed,  the  inertia  forces 
are  computed  and  added  to  the  applied  external  forces,  the  work  function 
is  formulated  and  tested  for  a  relative  minimum.  With  the  aid  of  the 
information  gained  in  this  teat,  a  new  configuration  is  found,  and  the 
process  is  repeated  until  the  equilibrium  imbalance  at  the  end  of  the 
time  step  is  sufficiently  small. 

This  solution  process  can  also  be  employed  to  obtain  the  nonlinear 
response  to  static  loads.  Aside  from  the  inertia  forces,  the  difference 
between  the  static  and  dynamic  analysis  is  conceptual.  Instead  of  a 
time  step,  a  load  increment  is  specified  and  the  corresponding  config¬ 
uration  is  again  obtained  by  work  function  minimization. 
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SECTION  2 


MATHEMATICAL  MODELS 

This  section  presents  mathematical  models  of  plane,  skeletal,  rein¬ 
forced  concrete  structures  and  their  environments. 

The  model  of  the  structure,  the  system  model,  is  a  discrete  model 
composed  of  line  elements  (models  of  reinforced  concrete  beam- columns) 
and  springs  (models  of  partial  joint  releases).  The  line  elements  admit 
geometric  and  physical  nonlinearities;  they  can  predict  the  behavior  of 
reinforced  concrete  members  subject  to  flexural  and  axial  distortions  up 
to  failure,  which  is  defined  as  the  limit  of  continuous  change  of  state. 

The  state  of  the  element  is  characterized  by  its  internal  energy.  The 
springs  are  restricted  to  linear  behavior. 

The  state  of  the  system  is  defined  by  the  work  function,  a  scalar 
function  that  contains  implicitly  all  the  forces  acting  on  the  system. 

The  work,  function  is  uniquely  defined  in  terms  of  the  generalized  co¬ 
ordinates,  which  must  be  related  to  the  equilibrium  path  (motion)  when 
the  system  behaves  nonlinearly  (cf.  section  2..2.4). 

Failure  criteria  are  formulated;  they  define  the  domain  in  which 
the  models  are  valid  and  provide  the  basis  for  predicting  element  fail¬ 
ure  and  structural  collapse. 

2.1  ACTIONS 

Actions,  mathematical  models  of  the  environment,  consist  of  the  self¬ 
weight  of  the  structure,  distributed  and  concentrated  loads,  inertia  forces, 
and  support  motions • 
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All  distributed  loads  and  9elf-weights  are  replaced  by  "equivalent" 
nodal  forces  [17].  In  the  linear  range  of  the  element,  the  equivalent 
nodal  forces  caused  by  transverse  member  loads  are  equal  in  magnitude 
and  opposite  In  sense  to  fixed-end  forces;  this  Is  a  consequence  of  the 
assumed  shape  functions  (cf.  section  2.2.1),  which  correspond  to  the 
homogeneous  solution  of  the  differential  equation  of  a  beam  In  flexure. 

This  property  does  not  exist  in  the  nonlinear  range  where  the  discrete 
element  forms  an  approximate  representation  of  the  continuum. 

Inertia  forces  are  computed  on  the  basis  of  lumped  masses  assigned 
to  the  nodal  degrees-of-freedom.  The  computation  of  the  lumped  masses 
follows  the  approach  described  in  reference  12  . 

2 . 2  ELEMENT  MODEL 

The  reinforced  concrete  beam-column  is  represented  by  a  gross  ele¬ 
ment  model.  This  means  that  the  element  forms  a  one-dimenBional  contin¬ 
uum,  which  is  discretized  in  the  modeling  process. 

The  initial  state  of  the  element  is  assumed  lo  be  unstrained.  Defor¬ 
mations  are  governed  by  the  fundamental  assumption  that  plane  sections 
remain  plane  and  normal  to  the  deformed  reference  axis.  Consequently, 
the  state  at  any  point  of  the  element  is  defined  by  the  state  of  the 
reference  axis.  Deformations  are  limited  by  the  assumption  that  strains 
and  rotations  are  small  relative  to  unity.  Axial  and  flexural  deformations 
are  modeled  explicitly;  only  a  measure  of  shear  distortions  and  their 
significance  is  provided.  Inelastic  deformations  are  modeled  up  to 
element  failure.  Structural  damping  is  incorporated  through  energy  dis¬ 
sipation  associated  with  inelastic  behavior. 

The  beam-column  effect,  the  coupling  of  axial  and  flexural  distortions. 
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is  represented  by  the  corresponding  nonlinear  term  in  the  strain-dis- 
placemenc  relation.  The  member-force  interactions ,  which  are  charac¬ 
terized  in  the  concrete  literature  by  behavio*.  models,  are  also  formu¬ 
lated  at  the  micro  level.  This  is  natural  since  the  behavior  model, 
a  macro  model  governing  the  axial  load -moment-curvature  relations  at  a 
section,  is  completely  defined  by  the  following  section  properties: 
the  strain  state,  the  constituents  of  the  section,  and  the  corresponding 
constitutive  laws.  The  variability  of  the  neutral  axis,  a  characteristic 
of  reinforced  concrete  beam-;  subjected  to  axial  and  flexural  distortions, 
is  modeled  by  admitting  axial  strain  variations  along  the  reference  axis. 
This  feature  is  illustrated  in  section  2.2.1. 

The  state  of  the  element  is  characterized  by  its  internal  energy. 
Conditions  of  equilibrium  are  formulated  for  the  assemblage  of  elements, 
the  structural  system.  The  modeling  process,  passing  from  the  continuum 
to  the  internal  energy  expressed  in  terms  of  a  finite  number  of  distor¬ 
tion  components,  is  depicted  schematically  in  Figure  3  :  u  and  v  define 

the  deformed  reference  axis;  u  is  a  4-dimensional  element  distortion 
vector  whose  components  represent  the  relative  element  displacements;  x 
and  y  are  the  coordinates  of  a  point  in  the  element  (Figure  5  ) ;  e  and 

a  denote  strain  and  stress  at  a  point,  respectively;  and  U  signifies  the 
internal  energy  of  the  element. 

2.2.1  DISCRETIZATION 

The  reference  axis  of  the  element  is  depicted  in  Figure  4  .  The 
reference  axis  must  lie  in  the  longitudinal  plane  of  symmetry,  the  plane 
of  bending,  of  the  element,  and  all  reference  axes  incident  at  a  joint 
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Fig.  3  MODELING  PROCESS 
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FIG.  4.  ELEMENT  REFERENCE  AXIS 


(a)  LONGITUDINAL  SECTION  (b)  CROSS  SECTION 

FIG.  5.  ELEMENT  SECTIONS 
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must  be  concurrent.  This  eliminates  the  modeling  of  joint  eccentricities. 
Moreover,  the  reference  axis  need  only  be  parallel  to  a  longitudinal  edge 
of  the  beam;  its  location  in  the  longitudinal  plane  of  symmetry  1b  arbi¬ 
trary  (see  illustrative  example  on  page  20). 

The  deformation  coordinate  axes,  the  x,  y-axes  in  Figure  4  ,  are 

defined  in  section  2.3.1.  The  transformation  of  the  global  joint  dis¬ 
placements  into  the  element  distortion  components  u^ ,  u^,  u^  is  presented 
in  section  2.3.1.  The  internal  distortion  compcnent,  u^, is  prescribed 
directly  in  the  solution  process. 

The  configuration  of  the  deformed  reference  axis  is  expressed  in  the 


form 

u(0  =  ^(Ouj  +  4>4(Ou4  (2.1) 

v(£)  -  <t»2(C)u2  +  ♦  (053  (2.2) 

where 

^  -  2C2  -  K  (2.3) 

*2  — 253  +  3£2  (2.4) 

*3  -  LU3  -  £2)  (2.5) 

\  -  4(-£2  +  O  (2.6) 


and 

S  *  x/L  (2.7) 

For  a  linear  element  satisfying  the  conditions  of  the  elementary 
flexure  theory,  Equation  2.2  represents  an  exact  description  of  the 
transverse  flexural  deflection  v  in  terms  of  the  relative  end-displace¬ 
ments  u2>  U3.  For  a  nonlinear  element,  the  shape  functions  <f>2,  $3 
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provide  only  an  approximate  representation  of  the  flexural  response. 

The  introduction  of  the  internal  distortion  component  u^  in  the  longi¬ 
tudinal  displacement  function, Equation  2.1,  permits  linear  variation  in 
the  normal  strain  along  the  reference  axis  (see  Equation  2.11).  This 
feature  makes  it  possible  to  describe  the  strain  state  corresponding  to 
a  linearly  varying  neutral  axis  with  respect  to  anv  reference  axis  in 
che  longitudinal  plane  of  symmetry.  This  property  is  illustrated  in  the 
following  example. 

Consider  the.  strain  state 
2e 

e(x,y)  =  -e - —  y  (i  +  — )  (2.8) 

°  h  L 

Oi  the  beam  shown  in  Figure  o  .  The  first  term  on  the  right-hand  side 
ot  Equation  2.8  represents  a  constant  normal  strain  induced  by  axial 
compression,  and  the  second  term  describes  a  flexural  strain  that  varies 
linearly  with  respect  to  the  orthogonal  reference  axes,  x  and  y;  the 
x-axis  coincides  with  the  centroidal  axis  of  the  beam;  h  and  L  denote 
the  height  and  length  of  the  beam,  respectively.  The  neutral  axis  is 
formed  by  the  straight  line  passing  through  points  P  and  Q  (Figure  6  ). 
introduce  a  reference  axis  that  does  not  coincide  with  the  centroidal 
axis;  e.g.,  let  the  location  of  the  reference  axis  be  described  by  the 
coordinate  transformations 


which  places  the  reference  axis  a  distance  h/4  below  the  centroidal 
axis  (Figure  6  ).  Substitution  of  Equation  2.9  into  Equation  2.8 
yields 
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CENTROIDAL  AXIS 


(2.10) 


‘-(x.  y)  =  -  lo  (j  -  ~)  -  2  _o  y(i  +  *) 

h 

A  comparison  oi  Equations  2.8  and  2.10  indicates  that  the  translation  of 
the  reference  axis  causes  the  normal  Btrain  to  vary  linearly  along  the 
reference  axis  but  does  not  alter  the  form  of  the  flexural  strain  term. 

Hence,  a  strain  state  corresponding  to  a  linearly  varying  neutral 
i:;is  can  he  described  rel ative.  to  areference  axis  that  admits  linearly 
varying  normal  strains. 


2.2.2  COMPATIBILITY 

The  point-wise  deformations  of  the  element  are  defined  by  the  strain- 


dis(  lacemcnt  relation  (Figure  3  ) 

2 

,  ,  du  1  .dv.  2  d  v 

<*,y>  -^+2  -y~2 

dx 


(2.11) 


where  e(x,v)  is  the  normal  strain  (in  the  x-direction)  at  any  point  (x,y) ; 
tl'.e  x-coordinate  locates  planes  nori^al  to  the  undeformed  reference  axis, 
and  the  y-coordinate  locates  points  in  that  plane;  u(x)  and  v(x)  define 
the  deflections  of  any  point  (x,0)  on  the  reference  axis  in  the  x  and  y 
directions,  respectively. 

The  terms  on  the  right-hand  side  of  Equation  2.11  admit  the  following 
geometric  interpretations:  The  first  term  defines  the  normal  strain 
induced  by  axial  deformations  of  the  reference  axis;  the  second  term 
represents  the  contribution  of  bending  of  the  reference  axis  to  the 
normal  strain  { 8 ] ;  i.e.,it  accounts  for  the  coupling  of  axial  and  flex¬ 
ural  distortions;  and  the  third  term  represents  the  elementary  bending 
s  train . 


J 
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Equation  2.11  is  valid  if  the  strains  and  rotations  are  small 
compared  to  unity  [e.g.,  3,  8,  11].  These  limitations  are  characteristic 
of  classical  stability  investigations  leading  to  conditions  of  infini¬ 


tesimal  stability  (e.g., the  Euler  buckling  load).  Equation  2.11  can 
form  the  basis  of  post-buckling  investigations  provided  the  strains 
remain  small  and  the  rotations  are  held  small  by  the  division  of  the 
element  into  sub-elements.  The  same  procedure  can  be  employed  to  model 
regions  of  large  distortions  induced  by  inelastic  deformations. 

With  the  aid  of  Equations  2.1  and  2.2,  the  normal  strain  can  be 
expressed  in  terms  of  the  element  distortion  components: 


U1  u4  1  uo  $3  _ 

c(5,n)  -  *;-  +  2  +ru3> 


3  -  \  2 


u2  *3 

•'“n’rV 


where 


(2.12) 


H  -  1 

(2.13) 

=  4(-2C  +  1) 

(2.14) 

*2  =  6(-C2  +  0 

(2.15) 

-  L(3£2  -  20 

(2.16) 

4>2  -  6(-2g  +  1) 

(2.17) 

^  =  L(6C  -  2) 

(2. 1C' 

and 


n  *  y/L 


(2.19) 
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2.2.3  CONSTITUTIVE  LAWS 


The  stress-strain  laws  governing  material  behavior  are  presented  in 
appendix  A.  They  are  expressed  in  terms  of  piece-wise  linear  functions 
sun  that  to  every  point  in  the  domain  (c)  corresponds  a  unique  point  in 
the  range  (a),  which  is  determined  by  the  strain  history. 

The  constitutive  laws  presented  model  the  behavior  of  concrete  (un- 
cont'ined  and  confined)  and  reinforcing  steel  for  monotonic  and  cyclic 
loading.  The  inherent  assumptions  and  limitations  are  stated. 

2.2.4  INTERNAL  ENERGY 

Energy  evaluation  represents  the  pivotal  task  in  the  search  of  the 
equilibrium  state  corresponding  to  a  prescribed  time  (or  load).  All 
measures  of  response  (e. g. .displacements,  deformations,  strains,  stresses, 
energies)  are  expressed  relative  to  the  "initial  state,"  which  is  the 
unstrained  and  unloaded  configuration  of  the  system. 

Energy  evaluation  in  the  context  of  the  solution  process  means  the 
computation  of  the  total  energy  of  the  system  for  a  given  displacement 
state.  The  internal  energy  evaluation  proceeds  as  follows:  On  the  basis 
of  Equation  2.11  and  appropriate  constitutive  laws  (appendix  A),  the 
"internal-energy  density,"  the  internal  energy  per  unit  volume,  is  de¬ 
termined.  Integration  of  the  internal-energy  density  over  the  volume  of 
the  element  yields  the  internal  energy  of  the  element.  The  internal 
energy  of  the  system  is  equal  to  the  sum  of  the  internal  energies  ol  all 
elements  comprising  the  system  (if  the  system  contains  release  springs, 
their  strain  energies  must  be  added). 

The  principal  assumption  in  the  energy  computation  is  that  no  "load 


reversals"  occur  during  a  time  step;  1. enduring  the  entire  time  step, 
the  strain  at  any  point  in  the  system  is  either  monotone  increasing  or 
monotone  decreasing. 

Energy  evaluations  must  be  conducted  numerically.  In  the  elastic 
range,  numerical  Integration  is  dictated  by  the  possible  variation  of 
cross-sectional  properties  (e.g.  area  of  compression  block)  over  some 
region  of  the  element.  For  instance,  an  axial  load  and  a  varying 
bending  moment  cause  a  varying  neutral  axis  (cf.  illustrative  example 
on  page  20  ) .  In  the  inelastic  range,  it  is  not  possible  to  formulate 
explicitly  the  variation  of  the  internal-energy  density  over  the  volume 
of  the  element. 

The  numerical  energy  evaluation  is  based  on  the  discretization  of 
the  energy  stored  in  the  element.  It  involves  two  principal  tasks: 

1.  The  computation  of  the  internal-energy  density  at  a  discrete 
number  of  points  in  the  element. 

2.  The  integration  of  the  internal-energy  density  over  the  volume 
of  the  element. 

The  computation  of  the  internal-energy  density  during  the  solution 
process  of  a  typical  time  step,  from  t^  to  is  described  with  the 
aid  of  Figure  7  .  t^  corresponds  to  the  time  at  which  the  last  equi¬ 
librium  state  of  the  system  has  been  obtained,  and  t£  denotes  the  tine 
at  which  the  next  equilibrium  state  is  sought.  The  stress-strain 
curves  in  Figure  7  govern  the  behavior  of  a  discrete  point  of  the 
element,  and  denote  strains  at  t^  and  t^»  respectively;  both 
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)  cases  are  illustrated.  The  inter- 


luniing  (t  )  and  unloading  (  V  *  . 

nal  •■iicd’v  densL’y  at  time  t0  is 

"*  ••  U*  d  H*  (2.20) 

»-  l  1 c. 


vhei 


(2.21) 


represents  the  internal-energy  density  at  t^,  and 


,  Jc 
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if 

if 


(2.22) 


represents  the  change  in  the  internal-energy  density  during  the  time 
step  tj,  ty.  It  follows  from  Figure  7  that  for  a  given  value  of  strain 
t, ,  there  corresponds  a  unique  value  of  stress.  Consequently,  the 
internal-energy  density. and  hence  the  internal  energy,  is  uniquely  defined 
bv  the  strain  state,  which  in  turn  is  a  unique  function  of  the  displace¬ 
ment  state.  Hence,  in  the  neighborhood  of  an  equilibrium  state,  the 
internal  energy  of  the  system  is  a  unique  function  of  the  generalized 
coordinates . 

In  the  inelastic  range,  the  internal-energy  density  consists  of  a 
dissipative  component  U*,  which  is  locked  into  the  material  by  residual 
stresses  on  the  microscopic  level,  and  a  recoverable  component  U*,  which 
is  released  by  the  ma.  ial  upon  unloading  (see  Figure  8  ).  The  dis¬ 
sipative  component  accounts  for  structural  damping. 
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(a)  LOADING 


(b)  UNLOADING 


FIG.  7.  INTERNAL- ENERGY  DENSITY 


(a)  LOADING  (b)  UNLOADING 

FIG.  8.  DISSIPATIVE  8  RECOVERABLE  ENERGY 
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The  computation  of  the  internal  energy  of  the  element  is  based  on 
the  Gaussian  quadrature  method  |  1 7 1 ;  the  concrete  and  steel  are  considered 
separately.  The  internal-energy  densities  are  evaluated  at  discrete 
points,  the  Gauss  points,  and  substituted  into  the  Gaussian  quadrature 
formula  to  yield  the  energy  stored  in  the  element.  The  Gauss  points  are 
distributed  in  the  longitudinal  plane  of  the  element  as  follows:  six 
points  (a  2x3  rule)  are  placed  in  the  ton  and  bottom  concrete  covers, 
nine  points  in  3x3  rule)  are  placed  in  the  concrete  between  the  covers, 
and  three  points  are  placed  along  the  centroidal  axis  of  each  steel  layer. 
The  accuracy  of  the  energy  computation  increases  with  the  number  of 
Gauss  points  per  element,  which  at  present  is  fixed.  Hence,  it  can  only 
be  controlled  indirectly  through  the  division  of  the  element  into  sub- 
elements  . 

t'nergv  variations  govern  the  behavior  of  the  mathematical  model  of 
the  structure.  The  accuracy  of  response  predictions  of  the  structure  is 
limited  by  the  accuracy  inherent  in  the  energy  evaluations.  For  this 
reason  internal  energies  induced  by  shear  distortions  are  not  included 
in  the  mathematical  model;  only  estimates  of  the  internal  energy  caused 
by  shear  distortions  and  measures  of  the  significance  of  these  distortions 
are  provided.  The  modification  of  the  element  model  to  account  for 
shear  distortions  introduces  uncertainties  which  may  seriously  affect  the 
reliability  of  the  model.  The  sources  of  uncertainty  are  identified 
in  the  following  discussion. 

On  the  basis  of  the  elementary  beam  theory  including  shear  effects, 
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the  shape  functions  In  liquation  2.2  can  he  modified  to  assume  the  I  orm  |I2| 


^2  -  t?7<-2«3  +  3f'2  +  YU 

and 


c2  +  \  a2  -e )] 


where 


12EI 

L2A  G 
s 


and 


(2.2  3) 


(2.24) 


(2.25) 


A  -  -  (2.26) 

S  K 

y  is  a  measure  of  the  relative  importance  of  shear  deformations.  In 
particular,  y  is  the  ratio  cf  the  shear  deflection  to  the  bending  de¬ 
flection  of  a  fixed-fixed  beam  subject  to  a  relative  end  displacement. 

It  is  important  to  recall  that  Equation  2.25  is  based  on  the  assumption 
that  the  beam  is  prismatic,  homogeneous,  isotropic,  and  linearly  elastic. 
Accordingly,  the  symbols  in  Equations  2.25  and  2.26  are  constants  for  a 
given  beam:  E  and  G  denote  Young’s  modulus  and  the  shear  modulus  cf 

elasticity,  respectively;  A,  A  ,  and  I  define  the  area,  the  effective 

s 

shear  area,  and  the  moment  of  inertia  of  the  cross  section,  respectively; 

k  is  a  shape  factor  that  reflects  the  variation  of  the  shear  stress 

across  the  section;  and  L  is  the  length  of  the  beam. 

For  an  inelastic  reinforced  concrete  beam-column,  the  quantities  in 

Equation  2.25  are  not  constants:  The  moduli  E  and  G  vary  pointwise  over 

the  volume  of  the  uncracked  concrete  and  steel;  the  section  properties 

A,  A  ,  and  I  vary  with  the  longitudinal  axis  of  the  beam  due  to  non- 
s 
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uniform  cracking;  especially  the  effective  shear  area  Ag  is  difficult 
to  define  since  the  shear-stress  distribution  over  a  cracked  section  is 
not  known.  In  essence,  the  problem  is  that  Equation  2.25  is  defined  in 
terns  of  macro  quantities  which  at  best  provide  an  indirect  description 
of  the  state  of  an  inelastic  reinforced  concrete  beam-column.  The  same 
difficulty  is  encountered  in  the  formulation  of  the  shear  energy  which 
is  defined  by  the  relation 


U 


T  2 
.  V~dx 

b  2A  G 
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(2. 27) 


or 


U 
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V21. 
2  A  G 


(2.28) 


s 

since  the  shear  force  V  is  constant  in  the  element  model. 

In  view  of  the  uncertainties  inherent  in  the  prediction  of  shear 
effects,  they  are  not  modeled  explicitly;  only  a  measure  of  the  sig¬ 
nificance  of  shear  distortions  is  provided  on  the  basis  of  Equation  2.25, 
and  an  estimate  of  the  internal  energy  induced  by  shear  distortions  is 
made  on  the  basis  of  Equation  2.28.  In  the  evaluation  of  Equations  2.25 
28,  E  and  G  are  assumed  to  be  elastic  and  ic  is  set  equal  to  1.20. 

2.2.5  STRESS  RESULTANTS 


* 


The  element  end-forces,  which  act  at  the  reference  axis  (cf.  Figure 
9  ),  are  computed  on  the  basis  of  the  following  formulas: 


f.  ,  *  f  o(L,y)dA 

A(L) 


(2.29) 
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-  ‘  yo(L,v)dA 
A(L) 


(2.30) 
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b3 


f  _  =  /  ycr(0,y)dA 

a3  A(0) 


(2.31) 


Oi.  b3  ai  bl  2 


fal  "fbl 


a2  b2 


where  fai’  fbi’  1  =  l>  2'  3’ 
respectively;  and  A(0) ,  A(L) 

a  &  b-end,  respectively. 


)/L  (2.32) 

(2.33) 

(2.34) 

are  the  element  forces  at  the  a  &  b-end, 
represent  the  cross-sectional  areas  at  the 


2.3  SYSTEM  MODEL 

The  system  model  is  a  mathematical  representation  of  plane,  skeletal, 
reinforced  concrete  structures.  It  is  an  assemblage  of  line  elements 
interconnected  at  a  finite  number  of  nodes.  The  elements  are  assumed  to 
be  rigidly  connected  at  the  nodes  unless  partial  or  complete  releases 
are  specified. 

In  the  linear  domain,  the  state  of  the  system  is  completely  defined 
in  terms  of  the  generalized  coordinates  which  consist  of  nodal  displace¬ 
ments,  internal-element  distortion  components,  and  relative  displacements 
at  releases.  In  the  nonlinear  domain,  the  generalized  coordinates  must 
ho  related  to  the  equilibrium  path  (motion)  of  the  system  to  define  the 
state  of  the  system  (see  section  2.2.4).  In  the  "initial  state,"  the 
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generalized  coordinates  are  zero. 

There  is  no  restriction  on  the  magnitude  of  the  generalized  coordi¬ 
nates  per  se;  however,  relative  displacements,  such  as  the  relative  dis¬ 
placements  of  nodes  linked  by  an  element,  are  limited  by  the  small  defor¬ 
mation  requirements  of  the  element  (cf.  section  2.2.2).  Violations  of 
these  limitations  can  be  resolved  through  the  insertion  of  additional 
nodes,  i.e., through  the  subdivision  of  elements. 

The  following  sections  are  concerned  with  compatibility  and  stability 
of  equilibrium  of  the  system. 

2.3.1  COMPATIBILITY 

This  section  relates  nodal  displacements  with  relative  element  dis¬ 
placements,  called  element  distortion  components.  In  the  derivation  of 
these  components,  four  orthogonal,  right-handed,  Cartesian  coordinate 
systems  are  employed;  they  are  called  global,  local,  joint,  and  deformation 
systems.  The  deflections  are  positive  if  they  take  place  in  the  positive 
direction  of  the  1,  2- axes;  the  positive  sense  of  rotations  about  the 
3-axis  is  determined  by  the  right-hand  rule. 

The  global  and  local  systems  correspond  to  the  coordinate  systems 
used  in  linear  matrix  analysis  (Figure  10).  Joint  coordinates  and 
joint  properties  (e.g., forces  and  displacements)  are  expressed  in  global 
coordinates  and  denoted  by  capital  letters.  Local  axes  define  the 
orientation  of  the  undeformed  element:  the  1-axis  coincides  with  the 
reference  axis,  and  the  2  &  3-axes  correspond  to  principal  axes  of  the 
cross-section.  The  1-axis  specifies  the  direction  of  the  element;  the 
element  in  Figure  10  goes  from  joint  i  to  joint  j.  Local  vectors  are 
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identified  by  lower-case  letters. 

The  transformation  of  a  two-dimensional  global  vector  Y  into  a 
two-dimensional  local  vector  y  is  defined  by  the  matrix  A: 


where 


and 


y  =  AY 


s 

c 


1 


c  =  cos  a,  s  =  sin  a 
It  follows  from  Fig.  10  that 
c  =  AX^/L ,  s  =  AX2/L 

where 


AX 


1 


(2.35) 

(2.36) 

(2.37) 

(2.38) 

(2.39) 


and  the  initial  element  length 
L  =  (AXX2  +  AX22)2 


(2.40) 


The  joint  and  deformation  reference  frames  are  moving  frames  of 
reference  rigidly  attached  to  the  joint  at  the  origin  of  the  element 
(Fig.  11).  In  the  initial  state,  the  joint  coordinate  system  coincides 
with  a  global  coordinate  system  originating  from  that  joint,  and  the 
deformation  coordinate  system  coincides  with  the  local  coordinate 
system. 

Vectors  expressed  in  joint  and  deformation  coordinate  systems 
are  identified  by  barred  capital  and  barred  lower-case  letters, 
respectively.  Since  the  joint  and  deformation  reference  frames  are 
fixed  relative  to  each  other,  corresponding  vectors  are  transformed 
by  the  matrix;  i.e.. 
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y  -  AY  (2.41) 

where  y  and  Y  are  vectors  expressed  In  deformation  and  joint  coordinates, 
respectively. 

The  global-joint  transformation  is  given  by 

Y  ■  BY  (2.42) 

where 

B  -[  ^  8ll  (2.43) 

"Si  Ci 

c^  »  cosU^  ,  *  sinU^  (2.44) 


and  is  the  rotation  of  joint  i  about  the  3-global  axis. 

If  follows  from  Eqs.  2.41  &  42  that  the  global-deformation 
transformation  is  defined  by 


y 


CY 


(2.45) 


where 


c  ■  ^  (2.46) 

The  derivation  of  the  element  distortion  components  follows 

directly  from  Fig.  11  .  The  relative  member-end  rotation 


U3  “  Uj3  -  Ui3 


(2.47) 


where  and  are  the  rotations  of  the  joints  j  and 
The  relative  member-end  deflections  u^,  are  expressed 


CAX*  -  d 


i,  respectively, 
in  matrix  form 
(2.48) 
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where 


u  -  [-1]  (2.49) 


AX* 

-  AX  + 

ah 

(2,0) 

AX 

■  vxi 

,  r  Xjl"Xil  , 

Xj2-Xi2  1 

AX 

■  '  ax'  ' 

(2.51) 

AU 

■  Y'l 

-  .  . 
°32-U12 

AU 

(2.52) 

and 

d  -  [  q  ]  =  AAX  (2.53) 

In  Eqs.  2.51-53,  X^,  X^  are  the  joint  position  vectors;  are  the 

joint  deflection  vectors:  and  d  defines  the  rigid-body  motion  of  the 
element.  With  the  aid  of  Eqs.  2.46,  50,  &  53,  Eq.  2.48  can  be 
reduced  to  a  form  suitable  for  numerical  evaluation: 


u  =  A(DAX  +  BAU) 


where 


n  =  B-I 


-s 


i 


s 


-2s 


(2.54) 


(2.55) 


si2  =  sin(Ui3/2) 


(2.56) 
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.nut  1  is  i.tu  Identity  matrix. 


For  infinitesimal  <1  i splncuments  (i.e., 


s . 


"tv  ci 


1 ,  s  |  j 


"nA", 


U.^  0),  To.  ?. .14  reduces  to 


u  =  A (FAX  +  U) 


(2.57' 


wlic  re 


i:  “  I 


n. 


n 


(2.58) 


2.3.2  STABILITY  OF  EQUILIBRIUM 

As  described  in  section  3.2,  the  search  for  the  equilibrium  state 
(  corresponding  to  a  set  of  prescribed  forces  is  governed  by  the  princinle 
o;  least  action;  i.e.,at  an  equilibrium  state  the  energy  function 
assumes  a  relative  minimum.  Thus,  if  an  equilibrium  state  is  found, 
it  is  a  stable  equilibrium  state. 

2.4  FAILURE  CRITERIA 

An  assemblage  of  elements  may  experience  element  and  system 
failure.  Fracture,  the  limit  of  continuous  change  of  state  [4], 
defines  element  failure.  System  failure  means  collapse  of  the 
assemblage . 

2.4.1  ELEMENT  FAILURE 

''Structure-sensitive"  properties  of  a  material,  such  as  the 
fracture  strength,  are  essentially  determined  by  local  imperfections 
in  the  groun  structure  of  the  material;  consequently,  they  exhibit 
•i  considerably  greater  degree  of  variability  than  "structure-insensitive" 
properties,  such  as  elastic  constants  [4],  Freudenthal  based  this 
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explanation  of  material  behavior  on  statistical  principles. 

Although  the  literature  reveals  significant  variations  in  the 
fracture  strength  of  concrete  and  reinforced  concrete  elements,  the 
corresponding  strength  criteria  are  seldom  based  on  probabilistic 
models;  i.e.,they  do  not  deal  with  these  inherent  uncertainties 
explicitly.  In  conventional  design,  the  problem  of  uncertain  failure 
strengths  is  usually  resolved  by  avoiding  such  failures  rather  than 
by  predicting  them.  The  underlying  philosophy  is  to  produce  ductile 
structures.  For  instance,  the  ultimate  moment  of  an  underreinforced 
concrete  beam  is  governed  by  the  yield  strength  of  the  steel.  Con¬ 
sequently,  the  significant  variability  of  the  crushing  strength  of 
the  concrete  has  little  affect  on  the  ultimate  flexural  strength  of  the 
reinforced  concrete  beam. 

In  this  project,  the  complete  structural  response  to  actions 
(including  system  failure)  must  be  predicted.  Under  static  actions, 
system  instability  without  element  failure  is  possible  (e.g.,the 
formation  of  a  collapse  mechanism)  but  perhaps  not  probable.  In 
the  dynamic  state,  it  may  not  be  possible  to  predict  the  collapse 
of  the  system  until  the  collapse  process  has  been  initiated,  in  which 
case  element  failure  is  probable.  In  any  event,  element  failure 
criteria  are  required. 

Since  fracture  appears  to  be  a  probabilistic  phenomenon  which 
is  not  modeled  explicitly,  it  is  monitored  via  lower-bound  criteria. 
When  the  possibility  of  fracture  is  detected,  the  user  must  decide 
whether  to  base  element  failure  on  the  conservative  lower-bound 
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criterion  or  to  modify  the  criterion  to  yield  more  probable  failure 
predictions  (see  appendix  B) .  This  procedure  requires  the  user  to 
recognize  and  deal  with  the  uncertainties  inherent  in  failure  criteria. 

Element  failure  criteria  are  resolved,  according  to  the  failure 
mechanisms,  into  micro  and  macro  criteria.  Micro  criteria  are  formulated 
on  the  basis  of  explicit  states  at  a  point,  such  as  the  strain  state. 
Macro  criteria  are  expressed  in  the  form  of  empirical  relations, 
involving  stress-resultants  and  element  properties. 

Micro  criteria  predict  primary  failures,  such  as  crushing  and 
cracking  of  concrete  and  fracture  of  steel,  induced  by  excessive 
normal  strains.  The  normal  strains  are  caused  by  flexural  and  axial 
distortions.  Crushing  of  concrete  may  occur  in  the  compression  zone 
of  unconfined  concrete;  it  may  also  take  place  in  conjunction  with 
compression  steel  "buckling"  in  confined  concrete.  Cracking  may  lead 
to  failure  if  it  initiates  in  an  unreinforced  region  of  a  beam  in 
flexure  or  if  the  entire  cross-section  is  in  tension.  Fracture  of 
steel  is  mainly  associated  with  very  light  reinforcement. 

Macro  criteria  are  concerned  with  shear-flexure  failures  [1] 
which  are  precipitated  by  the  formation  of  a  diagonal  tension  crack; 
the  resulting  failures  are  called  diagonal-tension,  shear-compression, 
and  shear-tension  failures.  The  nominal  average  shear  stress  is  used 
as  a  measure  of  the  diagonal  tension  strength.  For  unreinforced  webs 
the  occurrence  of  a  diagonal  tension  crack  is  regarded  as  element 
failure.  Although  diagonal  tension  cracks  tend  to  stabilize  in  short 
and  intermediate-length  beams,  the  crack  stabilization  mechanism  is 
not  well  enough  understood  to  warrant  utilization  of  the  reserve 
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strength  associated  with  shear-compression  and  shear-tension  failures. 
For  beams  with  appropriate  web  reinforcement,  the  web  reinforcement 
assures  the  stab i l 1 zat ion  of  the  diagonal  tension  crack;  however, 
yielding  of  the  web  reinforcement  can  lead  again  to  the  type  of  shear 
failures  experienced  by  the  unreinforced  beam. 

A  classification  of  all  possible  failure  modes  is  presented  in 
appendix  H.  In  addition,  lower  bound  criteria  are  stated,  and 
modifications  are  formulated  for  the  selection  of  more  probable 
failure  criteria. 

2.4.2  SYSTEM  FAILURE 

System  failure  can  be  linked  to  instability  of  equilibrium. 
Stability  is  the  property  of  equilibrium  to  sustain  disturbances. 

This  means  that  a  stable  system  remains  functional  in  the  perturbed 
state.  Degree  of  stability  of  equilibrium  is  a  measure  of  the  dis¬ 
turbances  an  equilibrium  state  can  sustain  [6]  .  If  an  equilibrium 
state  is  unstable  relative  to  a  particular  disturbance,  the  degree 
of  stability  is  zero. 

The  solution  process  employed  in  this  analysis  converges  only 
to  stable  equilibrium  states.  Hence,  the  problem  is  not  to  ascertain 
stability  of  equilibrium  but  to  predict  whether  an  equilibrium  state 
exists  for  a  prescribed  set  of  actions.  The  concept  of  degree  of 
stability  of  equilibrium  serves  as  a  basis  for  this  prediction.  The 
"average  curvature"  of  the  work  function  at  the  equilibrium  state  is 
selected  as  a  measure  of  degree  of  stability  of  equilibrium.  The 
computation  of  the  average  curvature  is  based  cn  the  values  of  the 
curvatures  of  the  work  function  at  the  equilibrium  state  in  the 
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direction  of  the  generalized  coordinates.  The  average  curvature  is 
not  likely  to  bo  zero  at  an  unstable  equilibrium  state  since  equilibrium 
is  unstable  if  the  minimum  principal  curvature  is  zero.  However,  the 
rate  of  change  of  a  load  parameter  with  respect  to  the  average  curvature 
approaches  zero  at  an  unstable  equilibrium  state.  Hence,  this  rate  of 
change  is  an  indicator  of  the  imminence  of  instability. 

The  relation  between  degree  of  stability  of  equilibrium  and  load 
level  is  depicted  in  Figure  12  ;  for  the  single-degree-of-freedom 

system,  the  curvature  of  the  work  function  at  the  equilibrium  state 
does  approach  zero  at  the  limit  load,  p*.  The  continuous  curve  over 
trie  domain  fKx<'x*  represents  stable  equilibrium  states,  and  the  broken 
curve  over  the  domain  x^x*  represents  unstable  equilibrium  states.  The 
decrease  in  degree  of  stability  of  equilibrium  with  increasing  load  is 
illustrated  by  the  work-function  curves  corresponding  to  the  equilibrium 
states,  x^,  x^,  x* .  The  respective  curvatures  at  the  equilibrium  points 
decrease  monotoni cally  to  zero.  For  a  load  in  excess  of  the  limit  load, 
e.g.,p  =  p^,  no  equilibrium  state  exists,  and  the  solution  process  employed 
in  this  study  cannot  converge. 

As  the  unstable  equilibrium  state  of  the  system  is  approached, 
a  load  increment  could  easily  push  the  load  beyond  the  limit  load. 

To  prevert  a  lengthy  search  for  an  equilibrium  state  that  does  not 
exist,  the  solution  process  is  terminated  after  the  deviation  from 
the  last  equilibrium  state  exceeds  a  prescribed  bound. 


42 


43 


i.rn  i'ations 


Tin*  prim:  i pa J  IImit:ntlmiK  and  approx Imat  ! on  -t  of  t  he  mathematical 
moilel  of  the  skeletal  reinforced  concrete  structure  are  summarized 
below ; 


1.  The  element  model  is  subject  to  the  standard  limitations 
associated  with  the  discretization  approach  of  the  finite 
element  method  (e . p internal  element  displacements  are 
expressed  approximately  in  terms  of  the  nodal  displacements; 
distributed  loads  are  replaced  by  "equivalent"  nodal  forces). 

2.  Plane  sections  arc  assumed  to  remain  plane  and  normal  to  the 
deformed  reference  axis  of  the  reinforced  concrete  beam. 

This  assumption  appears  to  be  reasonable  up  to  the  formation 
of  diagonal  tension  cracks  of  unreinforced  webs  [16],  which 
represent  limits  of  continuous  change  of  state  of  the  element. 

3.  Normal  strains  and  rotations  are  assumed  to  be  "small"  in 
the  sense  that  their  squares  are  negligible  with  respect  to 
unity  [3,  111;  i.e.,  they  are  regarded  to  be  infinitesimals. 
These  limitations  are  acceptable  since  the  fracture  strains 
of  the  materials  modeled  in  this  project  meet  this  require¬ 
ment,  and  the  rotations  can  be  controlled  through  element 
subdivision.  Shear  distortions  are  not  modeled  explicitly; 
the  indication  is  that  a  modification  of  the  gross  element 
model  to  include  shear  deformations  is  likely  to  impair  the 
quality  of  the  model. 

4.  The  constitutive  laws  governing  material  behavior  are  described 
by  deterministic  models.  Consequently,  they  represent  at 
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best  Che  statistical  mean  of  the  material  properties  and 
do  not  reflect  the  significant  randomness  characteristic 
of  some  properties  such  as  the  fracture  strength. 

5.  Energy  computation  is  based  on  the  assumption  that  no  "load 
reversals"  occur  during  a  time  (or  load)  increment  of  the 
solution  process.  Moreover,  the  computation  of  the  internal 
energy  is  based  on  the  evaluation  of  the  internal  energy 
densities  at  a  discrete  number  of  points  in  the  beam  element. 
This  introduces  another  discretization  error,  which  vanishes 
only  in  the  limit. 

6.  Element  failures  precipitated  by  material  fractures  are 
inherently  random  phenomena  which  can  only  be  monitored  by 
lower  bound  criteria  in  a  deterministic  analysis.  In  the 
event  that  structural  collapse  is  strongly  influenced  by 
element  failures  (as  in  contrast  to  the  formation  of  a 
"plastic"  collapse  mechanism),  the  quality  of  this  prediction 
by  deterministic  methods  is  questionable. 
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MOTION  1 


RESPONSE 

The  response  of  the  system  model  to  act  ons  is  sought  at  a  discrete 
niinihor  of  points  in  time.  The  solution  process,  formulated  by  Melosh  and 
Kelley  (*»],  is  a  closed  iterative  process  within  two  successive  time 
points:  'rhc  state  of  the  system  is  assumed  to  be  known  at  the  beginning 
and  is  sought  at  the  end  of  the  time  step.  Thus,  if  the  state  of  the 
system  is  known  at  one  point  ir  time,  the  response  determination 
proceeds  like  a  chain  reaction  through  successive  discrete  points. 

The  solution  process  comprises  two  fundamental  concepts: 

1.  discretization  of  motion 
?.  work-function  minimization. 

The  motion,  time  functions  of  the.  generalized  system  coordinates,  is 
discretized  via  the  finite  element  method  ( 17 ] ;  this  process  is  analogous 
to  N’ewmark's  : -method  ( 10 ]  :  Each  displacement  function  is  completely 
defined  over  a  time  step  by  three  initial  conditions,  which  are  known, 
and  one  end  condition,  which  is  the  desired  displacement  at  the  end  of 
the  time  step.  The  work  function  [7],  a  scalar  function  that  contains 
implicitly  all  the  forces  acting  on  the  system  (applied,  inertia, 
internal)  ,  is  expressed  in  terms  of  the  unknown  system  coordinates  at 
the  end  of  the  time  step.  The  desired  system  configuration  is  obtained 
by  minimization  of  the  work  function,  which  assumes  a  relative  minimum 
at  the  dynamic  equilibrium  state. 

Function  minimization  is  based  on  Stewart's  modification  of 
Davidson's  method  [13].  A  measure  of  the  quality  of  the  response 
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predictions  Is  provided  through  error  controls  linked  to  automatic 
time-step  selections. 

3.1  DISCRETIZATION  OF  MOTION 

The  time  domain  Is  subdivided  Into  time  segments  At,  and  the  dis¬ 
placement  functions  are  approximated  over  each  subdomain  by  a  finite 
power  series  of  the  form 

X1  (t)  -  a1Q  +  a^t  +  «l2t2/2  +  a13t3/6,  0±t<At  (3.1) 

where  x^  represents  the  1^  generalized  system  coordinate,  and  t  is  the 
normalized  time  coordinate.  The  constant  coefficients  in  Eq.  3.1 
are  determined  on  the  basis  of  the  following  end  conditions: 


x  t 

ai 

-  xl(0) 

(3.2a) 

• 

x  * 
ai 

'  3txi(0) 

(3.2b) 

2 

Xai 

■  — »x  (0) 
dt  x 

(3.2c) 

*bi 

■  xt  (At) 

(3. 2d) 

where  x  .,  x  .,  x  .  denote  the  displacement,  velocity,  acceleration, 
ai  ai  ai 

respectively,  at  the  beginning  of  the  time  step,  and  x^  denotes  the 
displacement  at  the  end  of  the  time  step.  It  follows  from  Eqs.  3.1  and 
3.2  that  the  displacement  and  acceleration  functions  can  be  expressed 
over  the  domain  [0,  At]  in  the  form 

x  (t)  ■  x  +  x  t  +  x  ,t2/2  +  6  t3/6  (3.3) 

i  ai  si  ai  i 
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X  (t)  *=  X  +  fi  t 

i  ai  i 


(3. A) 


where 


*1  "  6(Kbi  * 


(3.5) 


j .  2  WORK-FUNCTION  MINIMIZATION 

Conditions  of  dynamic  equilibrium  are  established  on  the  basis  of 
the  principle  of  virtual  work,  which  states  that  the  vanishing  of  the 
virtual  work  for  all  possible  virtual  displacements  represents  a 
sufficient  condition  of  equilibrium;  i.e., 


6W  =  0 


(3.6) 


for  all  independent  virtual  displacements  is  a  sufficient  condition  of 
equilibrium.  The  total  virtual  work  can  be  expressed  as 


4W  =  6W  -  50 
e 


(3.7) 


where  5W^  represents  the  virtual  work  of  external  forces  and  6U  denotes 
the  first  variation  in  the  internal  energy  of  the  system: 


T 

5xb  P1 


(3.8) 


6U  =  5x,  r. 


(3.9) 


Pb  =  fb  +  fb  '  “b 


(3.10) 


and  the  i  component  of  r^  is 


bi  ^x^ (At) 


(3.11) 
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The  subscript  b  In  Gqs.  3.8-11  signifies  that  the  corresponding  variables 
are  evaluated  at  the  end  of  the  time  step,  at  t  ■  At;  ,  x^  denote 
the  generalized  displacement,  acceleration  vectors, respectively;  m  is 
a  diagonal  mass  matrix;  p.  is  the  generalized  external  force  vector, 

D 

g 

which  consists  of  the  applied,  f^,  equivalent,  f^,  and  Inertia,  -mx^, 
force  vectors;  r^  represents  the  generalized  Internal  force  vector, 
whose  components  are  the  partial  derivatives  of  the  internal  energy 
with  respect  to  the  generalized  coordinates.  The  superscript  T 
signifies  transposition.  Eqs.  3.6-9  lead  to  the  condition 


6xb  (pb  '  V  "  0 


(3.12) 


which  yields  the  equilibrium  equation 


Pb‘rb  "  0 


(3.13) 


In  the  vicinity  of  the  equilibrium  state  corresponding  to  the 
beginning  of  the  time  step,  U  is  a  function  of  the  generalized 
coordinates  (cf.  section  2.2.4).  Moreover  p^  is  a  function  of 
by  virtue  of  Eqs.  3.4,  5.  Thus,  the  equilibrium  equation,  Eq.  3.13, 
is  a  function  of  x^. 

The  unknown  generalized  coordinates  at  the  end  of  the  time  step 
are  not  obtained  by  direct  solution  of  Eq.  3.13  but  by  minimization 
of  the  corresponding  work  function 

W(xfe)  »  x£pb  -  iKx^  (3.14) 
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The  sut ionary  condition,  Eq .  3.6,  which  leads  to  the  equation  of  dynamic 
equilibrium  is  also  a  minimum  condition.  On  the  basis  of  the  principle 
of  least  action  [/],  the  work  function  W  assumes  a  relative  minimum 
at  x^,  t he  solution  of  Kq.  3.13. 

3.3  PkOCFSS  ERRORS 

There  are  essentially  two  sources  of  error  in  the  solution  process  [?]: 
truncation  error  and  iteration  error.  The  truncation  error  is  induced  by  the 
approximate  representation  of  the  displacement  function  over  a  time  step 
by  a  finite  power  series.  The  truncation  error  decreases  with  the  size 
of  the  time  step  and  vanishes  in  the  limit;  hence,  it  can  be  controlled 
by  varying  the  length  of  the  time  step.  The  iteration  error  arises  in 
the  minimization  process,  which  converges  in  the  limit  to  the  exact 
solution.  Hence,  the  iteration  error  can  be  made  arbitrarily  small 
by  a  sufficiently  large  number  of  iterations. 

The  force  imbalance  at  the  mid-point  of  the  time  step  is  selected 
as  a  basis  for  a  measure  of  the  truncation  error.  It  follows  from 
Eq.  3.13  that  the  unbalanced  iC^  generalized  force  component 

^(t)  “  Pi(t)  -  r^ (t)  ,  0<t<At  (3.15) 

The  relation 

e  (t)  =  [Mt)x ,(t)]/W(t)  (3.16) 

i  i  i 

transforms  the  force  imbalance  into  a  relative  energy  imbalance. 
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Denote 


efi  -  max | e^(0) j  (3.17a) 

eab  -  max|ei(At/2)|  (3.17b) 

■  max  |  ^(At)  | ,  i  -  1,  2,  ...»  n  (3.17c) 


ea>eb  constitute  measures  of  the  Iteration  error  at  the  beginning  and 
end  of  the  time  step,  respectively,  and  eab  is  a  measure  of  the  truncation 
and  iteration  errors  at  the  mld-polnt;  n  Is  the  number  of  generalized 
coordinates.  If  one  assumes  that  the  iteration  error  varies  linearly 
over  the  time  step,  a  measure  of  the  truncation  error  is  obtained  in 
the  form  (cf.  Fig.  13  ) 

*T  '  \b  -  <«.  +  %)/2  (3-18> 


The  length  of  the  time  step  is  governed  by  the  following  inequality 


e 


<  e<e 
a  —  t-  u 


(3.19) 


where  e „  and  e  define  lower  and  upper  bounds  on  the  truncation  error 
measure,  respectively.  The  time  step  is  increased  i.'  e  <e  and  decreased 

X  X 

if  eT>eu-  The  lower  bound  is  imposed  to  assure  computer  accuracy; 
i.e.,  to  assure  that  the  time  step  is  large  enough  to  produce  measurable 
changes  in  the  response.  The  relation  between  computer  error,  truncation 
error,  and  step  length  is  depicted  in  Fig.  14  .  The  accuracy  of  the 
solution  process  is  apparently  insensitive  to  variations  in  At  over 
the  domain  (At^,  At^) .  The  most  economical  step  is  near  At^> 
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FIG.  13  TRUNCATION  »  ITERATION  ERRORS 


5? 


4  < 


A  direct  measure  of  the  iteration  error  la  provided  by  the  maximum 
absolute  value  of  the  unbalanced  generalized  force  component  at  the  end 


of  the  time  step 

ll>b  ■  max|ipi(At)  |  ,  1  -  1,2 . n 

The  minimization  process  is  continued  until 


*b  i». 


where  tpu  Is  a  prescribed  upper  bound  on  the  force  Imbalance. 


(3.20) 
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SECTION  4 


SUMMARY 

This  report  describes  the  mathematical  models  and  the  solution 
process  which  form  the  basis  of  the  computer  program  SINGER.  The 
function  of  SINGER  is  to  predict  the  behavior  of  plane  skeletal  rein¬ 
forced  concrete  structures  in  their  environments.  Of  primary  interest 
is  the  transient  nonlinear  response  Including  element  failures  and 
structural  collapse. 

The  principal  features  of  the  mathematical  models  and  the  solution 
process  are  summarized  below: 

ACTIONS 

Actions,  mathematical  models  of  the  environment,  consist  of  the 
self-weight  of  the  structure,  distributed  and  concentrated  static  and 
dynamic  loads,  inertia  forces,  and  support  motions.  All  distributed 
forces  are  replaced  by  equivalent  nodal  forces.  Lumped  masses  are 
assigned  to  the  nodal  degrees-of-f reedom. 

SYSTEM  MOD7'! 

a.  The  structure  is  represented  by  an  assemblage  of  line  elements 
(models  of  reinforced  concrete  beam-columns)  and  springs  (models 
of  partial  joint  releases)  interconnected  at  a  finite  number  of 
nodes. 

b.  The  state  of  the  system  is  characterized  by  the  work  function, 

a  scalar  function  that  contains  implicitly  all  the  forces  acting 
on  the  system.  The  work  function  is  uniquely  defined  in  terms 
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of  Che  generalized  coordinates,  vhlch  must  be  related  to  the 
equilibrium  path  (motion)  whan  the  system  behaves  nonllnearly. 

c.  There  Is  no  direct  restriction  on  the  magnitude  of  the  generalized 
coordinates,  which  consist  of  nodal-displacements,  relative  release- 
displacements,  and  Internal  element-displacements.  However,  relative 
displacements  of  nodes  linked  by  elements  are  limited  by  the  small 
deformation  requirements  of  the  elements.  Violations  of  these 
limitations  can  be  resolved  through  subdivision  of  the  elements. 

d.  The  transformation  of  the  large  nodal  displacements  Into  relative 
element  displacements  is  expressed  with  the  aid  of  two  frames  of 
reference:  The  global  frame  of  reference,  which  is  fixed  In 
space,  Is  used  to  describe  nodal  properties  (e.g., initial  state, 
displacements,  forces);  the  deformation  frame  of  reference,  a 
moving  frame  of  reference,  is  used  to  describe  element  properties 
(e.g. .strains,  stresses,  distortions). 

e.  In  a  static  analysis,  system  failure,  structural  collapse,  is  linked 
to  Instability  of  equilibrium.  In  a  dynamic  analysis,  structural 
collapse  is  inferred  from  the  motion  of  the  system. 

ELEMENT  MODEL 

a.  The  beam-column,  the  basic  structural  element,  is  modeled  as  a  one¬ 
dimensional  continuum,  which  is  discretized.  Axial  and  flexural 
deformations  are  modeled  explicitly;  only  a  measure  of  shear 
distortions  and  their  significance  is  provided.  Deformations  are 
limited  by  the  assumption  that  strains  and  rotations  are  small 
relative  to  unity.  Inelastic  deformations  are  modeled  up  to  element 
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failure.  Energy  dissipation  induced  by  inelastic  behavior  accounts 
for  structural  damping. 

b.  The  beam-column  effect,  the  coupling  of  axial  and  flexural 
deformations  is  represented  by  the  corresponding  nonlinear  term 
in  the  strain-displacement  relation.  The  varying  neutral  axis, 
a  characteristic  of  beam-columns,  is  modeled  by  admitting  normal 
strain  variations  along  the  reference  axis.  This  feature  makes 
it  also  possible  to  locate  the  reference  axis  anywhere  in  the 
longitudinal  plane  of  symmetry  of  the  element;  thus  it  eliminates 
modeling  of  joint  eccentricities. 

c.  Excessive  deformations  associated  with  slender  elements  or  "plastic" 
hinges  are  controlled  by  the  division  of  the  element  into  subelements. 

d.  Constitutive  laws  for  concrete  (unconfined  and  confined)  and  rein¬ 
forcing  steel  are  described  in  the  form  of  piecewise  linear  stress- 
strain  curves.  Material  behavior  under  monotonic  and  cyclic 
loading  is  modeled. 

e.  Element  failure,  which  is  defined  as  the  limit  of  continuous  change 
of  state,  is  predicted  on  the  basis  of  lower-bound  criteria. 
Modifications  of  these  criteria  are  formulated  to  permit  more 
probable  failure  predictions. 

RESPONSE 

The  solution  process  initiates  at  a  point  where  the  state  of  the 
system  is  completely  defined  and  proceeds  along  discrete  points  of 
the  motion:  The  time  function  of  each  generalized  coordinate  is  approxi¬ 
mated  within  two  successive  points  in  time,  the  time  step,  by  a  finite 
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power  series  whose  undetermined  coefficient  corresponds  to  the  unknown 
displacement  at  the  end  of  the  time  step.  This  representation  of  the 
motion  permits  one  to  formulate  the  work  function  of  the  system  at  the 
end  of  the  time  step  in  terms  of  the  unknown  generalized  coordinates. 

The  desired  equilibrium  state  is  obtained  by  minimization  of  the  work 
function. 

LIMITATIONS 

Spatial  and  temporal  discretization  and  the  inherent  variability 
of  material  properties  form  the  principal  sources  of  error.  Spatial 
discretization  errors  can  be  controlled  through  the  subdivision  of 
elements.  Although  internal  energy  computation  is  based  on  a  fixed 
mesh  imposed  on  the  longitudinal  plane  of  the  element,  the  reduction 
of  the  element  length  results  in  a  mesh  refinement,  and,  hence  it 
improves  the  accuracy  of  the  energy  computation.  Temporal  discretization 
errors  can  be  controlled  by  varying  the  size  of  the  time  step.  The 
constitutive  laws  governing  material  behavior  are  described  by  determin¬ 
istic  models,  which  do  not  reflect  the  randomness  of  some  properties  such 
as  the  fracture  strength.  Consequently,  element  failures  precipitated 
by  material  fracture  are  monitored  via  lower-bound  criteria. 
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NOTATION 


A, A 
’  s 

A ,  B ,  C ,  D ,  F 
B,B 

c 

C.i 

d 


e 


i 


e 


a,eb 


'ah 


e,r 


e„  ,e 

i  u 

E,G 


f  .fu 
a  b 


£al’fbi 

Vfb 


h 

I 


*r 


area,  effective  shear  area  of  beam 
transformation  matrices 
compatibility  matrices 
cos  a 

cos 

rigid-body  configuration  of  element  in  deformation 
coordinates 

energy  imbalance  corresponding  to  generalized  force 
component 

measures  of  iteration  error  at  the  beginning,  end  of 
time  step  At 

measure  of  truncation  and  iteration  errors  at  At/2 
measure  of  truncation  error 
lower,  upper  bound  on  e^, 

Young' s;  shear  modulus  of  elasticity 
element  force  vectors  at  the  a,  b-end 
element  forces  at  the  a,  b-end 

generalized  applied,  equivalent  force  vector  at  the 

end  of  time  step  At 

height  of  beam 

moment  of  inertia  of  beam 

identity  matrix 

stiffness  coefficients 

tangent  stiffness  matrix 
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L  « 


length  of  bean 


diagonal  mass  matrix 

number  of  generalized  coordinates 

generalized  external  force  vector,  at  the  end  of  time 

step  At 

unbalanced  force  vector  (scalar) 

generalized  Internal  force  vector  at  the  end  of  time 
step  At 
sin  a 


specific  values  of  t 
time  atep 

deflections  of  point  (x,0)  on  the  reference  axis 
element  distortion  vector,  component 
internal  energy 


1st  variation  of  internal  energy 

internal  energy  induced  by  shear  deformations 

displacement  vector  of  joint  1  in  global  coordinates 

rotation  of  joint  1  about  3-global  axis 

relative  joint  displacement  vector  in  global  coordinates 

internal-energy  density 

dissipative,  recoverable  internal-energy  density 
internal-energy  density  at  time  t^,  t2 
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II* 

12 


V 

V 

w 
sw 
f  w 

e 

x,y 


x,x. 


lx 


V*b 


X  .  ,x  . ,x  . 
al  ai  al 


*bi 

6xb 

VXj 

XU 

AX=X^-X1 
AX1 ,AX2 
AX*=AX+AU 


ci 


Y 


change  in  inLernal-energy  density  during  the  time 
interval  t^, 
shear  force 
volume 

work  function 

total  virtual  work 

external  virtual  work 

element  deformation  axes 

generalized  coordinate  vector,  component 

nC^  trial  solution 

correction  to  x 

n 

generalized  displacement,  acceleration  vector  at  the  end  of 
time  step  At 

displacement,  velocity,  acceleration  at  the  beginning  of  time 
step  At 

displacement  at  the  end  of  time  step  At 
virtual  displacement  vector  at  the  end  of  time  step  At 
position  vectors  of  joints  i,  j  in  global  coordinates 
global  coordinate  of  joint  i  in  j  direction 
relative  position  vector  in  the  initial  state 
components  of  AX 

relative  position  vector  in  the  displaced  state 

angle  between  1-local  axis  of  element  and  1-global  axis 

coefficient 

shear  deformation  factor 
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0 


K 

a 


virtual  variation 

normal  strains 

normal  strain  at  time  t^ 

nondimens lonal  element  deformation  axes 

shape  factor 

normal  stress 

element  shape  function 

d4>A 

dT 

d24>i 


measure  of  iteration  error  at  the  end  of  time  step  At 

1th  unbalanced  generalized  force 

upper  bound  on  ik 
b 
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APPENDIX  A 


CONSTITUTIVE  LAWS  FOR  CONCRETE  AND  STEEL 

This  appendix  summarizes  the  uniaxial  stress-strain  curves  used  in 
describing  the  material  response  of  a  single  fiber  of  either  concrete  or 
steel.  The  curves  described  are  the  default  stress-strain  curves  gener¬ 
ated  bv  the  program.  The  user  has  the  option  to  specify  others  if  he  so 
desires . 

A. 1  ASSUMPTIONS  ON  MATERIAL  BEHAVIOR 

The  following  assumptions  have  been  made  in  developing  the  consti¬ 
tutive  models  presented  herein: 

1.  Stresses  in  the  concrete  and  steel  are  uniquely  related  to  the 
strains.  For  direct  tension  and  compression  tests  under  short 
time  loading,  this  is  correct.  This  permits  the  calculation  of 
the  stresses  in  the  concrete  and  steel  once  the  strains  are 
known . 

2.  The  stress-strain  relationship  for  compressed  concrete  not 
confined  by  lateral  reinforcement  is  identical  to  that  for 
concrete  in  direct  compression.  The  neglect  of  a  strain 
gradient  effect  in  the  compression  zone  of  a  beam  is 
justified  by  adequate  correlation  between  experimental 
results  and  many  flexural  theories  based  on  this  assumption. 

3.  The  stress-strain  relationship  for  compressed  concrete  confined 
by  lateral  reinforcement  has  a  strength  greater  than  the  un¬ 
confined  direct  compression  strength.  Data  are  presented  in 
section  A. 2  which  supports  and  quantifies  this  assumption. 
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4.  Tension  stress  In  concrete  Is  neglect**'?.  The  magnitude  of  the 
tensile  stresses  In  the  concrete  Is  small  compared  to  that  In 
the  reinforcement  and  their  neglect  will  not  significantly 
change  the  results  of  the  analysis. 

5.  Concrete  stress-strain  curves  are  valid  for  normal  weight 
concrete  and  compressive  strengths  between  2500  psl  and  8000 
psl.  Lightweight  and  heavyweight  concretes  are  excluded. 

There  Is  sufficient  test  data  to  generalize  the  curves  pre¬ 
sented  to  only  a  limited  range  of  concrete  strengths. 

6.  The  stress-strain  relationships  for  steel  can  be  determined 
from  tension  tests  for  both  the  behavior  In  tension  and 
compression.  Complete  stress-strain  curves,  Including  strain 
hardening  and  breaking  strengths,  are  given  in  section  A. 3. 

These  are  limited  to  steels  with  yield  points  from  33  ksl  to 
75  ksl. 

7.  Creep,  shrinkage,  and  temperature  effects  are  ignored.  For 
short  duration  loadings,  the  first  two  effects  can  be  neglected. 
The  change  of  material  properties  with  temperature  is  not 
sufficiently  documented  for  reinforced  concrete  and  therefore 

is  omitted. 

8.  Strain  rate  effects  on  material  response  are  neglected.  This 
can  influence  the  stress-strain  response  at  local  points  in 

the  structure.  However,  it  is  assumed  that  the  overall  response 
of  the  structure  will  not  be  significantly  affected  by  ignoring 
this  complexity. 

9.  Adequate  lateral  support  is  present  to  prevent  buckling  of 
steel  in  compression.  This  assumption  is  valid  as  long  as 
the  concrete  cover  is  Intact.  After  spalling  has  taken  place, 
lateral  support  must  be  provided  by  lateral  ties  or  stirrups. 
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In  developing  the  computer  code,  checks  are  made  on  the  above  assump¬ 
tions  whenever  possible.  For  example,  the  concrete  compressive  strength 
given  must  be  within  the  range  specified,  stirrup  spacing  is  checked 
against  a  requirement  for  prevention  of  local  buckling,  etc.  If  one  of 
the  assumptions  is  violated,  a  warning  is  given  to  the  user  that  he  is 
using  the  program  beyond  its  intended  application. 

A.  2  STRESS-STRAIN  RELATIONSHIP  FOR  CONCRETE 

If  concrete  is  compressed  in  one  direction,  it  tends  to  expand 
laterally.  If  this  expansion  occurs  freely,  the  concrete  is  said  to  be 
"unconf ined"  and  principal  compressive  stresses  exist  in  one  direction 
only.  On  the  other  hand,  if  such  lateral  expansion  is  restricted,  the 
concrete  is  said  to  be  "confined"  and,  as  a  result  of  such  restriction, 
compressive  stresses  develop  in  all  directions.  Up  to  the  stage  corre¬ 
sponding  to  crushing  the  behavior  of  the  concrete  is  essentially  that 
of  the  unconf  ined  concrete.  Beyond  this  stage,  the  concrete  core  bound 
by  lateral  reinforcement  has  greater  strength  and  ductility  than  the  un¬ 
confined  concrete.  Because  of  these  differences,  it  is  necessary  to 
describe  stress-strain  relationships  for  both  types  of  concrete. 

A. 2.1  Stress-Strain  Curve  for  Unconfined  Concrete 

The  default  stress-strain  curve  for  unconfined  concrete  is  given  in 

Figure  A.l.  A  non-dimensional  plot  is  not  possible  because  the  slope  of 

the  descending  branch  is  dependent  on  the  compressive  cylinder  strength, 

f'.  The  curve  is  divided  into  two  portions,  AB  and  BC'.  For  the  region 
c 

AB,  a  parabolic  expression  (represented  by  a  series  of  straight  line 
segments  in  the  program)  given  in  reference  A.l  is  used: 
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FIG.  A.  h  COMPLETE  STRESS- STRAIN  CURVE  FOR  CONCRETE 


(A.  1) 


in  which  the  strain  at  maximum  stress  is  assumed  to  be  e  ■  0.002.  It  la 

o 

also  assumed  that  the  maximum  stress  Is  the  cylinder  strength  f'  ,  l.e., 

c 

the  factor  0.85  Is  not  Included.  The  reason  for  this  is  that  the  0.85 

factor  was  based  on  column  tests  without  a  strain  gradient.  When  a 

strain  gradient  is  present,  such  as  in  a  member  in  bending,  observations 

have  shown  (reference  A. 2)  that  a  factor  of  1.0  is  conservative.  The  region 

BC'  is  defined  by  a  straight  line  whose  slope  is  determined  by  the  strain 

Sou’  w^en  c^e  concrete  stress  has  fallen  to  50%  of  the  cylinder  strength 

of  the  unconfined  concrete.  This  is  given  in  reference  A. 3  as 

3  +  0.002f ' 
c 

e50u  *  f'  -  1,000  (A*2) 

c  * 

in  which  is  expressed  in  pounds  per  square  inch.  The  straight  line  is 
continued  until  the  concrete  strain  reaches  the  failure  value  defined  in 
Appendix  B,  section  B.3.1,  At  this  point  the  unconfined  concrete  is  no 
longer  effective  and  is  removed  from  the  cross-section. 

A. 2. 2  Stress-Strain  Curve  for  Confined  Concrete 

The  default  stress-strain  curve  for  confined  concrete  is  also  given 
in  Figure  A.l.  The  curve  is  divided  into  four  regions  and  is  similar  to 
the  curve  given  in  reference  A. 4.  For  the  region  AB,  the  curve  is  identical 
to  that  given  by  Eq.  A.l. 

For  the  region  BB',  a  modification  in  the  curve  of  reference  A. 4  is  made 
to  include  an  increase  in  compressive  strength  when  lateral  ties  are 
present.  Recommended  increases  for  this  region  vary  from  nearly  50%  (reference 
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A. 5)  to  zero  (reference  A. 6).  However,  a  majority  of  the  researchers  indicate 
that  a  modest  increase  is  reasonable,  and  the  following  expression  is  used: 

f"  -  f'  +  Af  (A. 3) 

c  c  c 

in  which  f"  -  confined  concrete  compressive  strength  and  Af  -  increase  in 
c  c 

compressive  strength  over  the  unconfined  value.  The  magnitude  of  Af'  is 

c 

dependent  on  the  confining  action  of  the  transverse  reinforcement.  A  the¬ 
oretical  discussion  in  reference  A. 5  indicates  that  the  lateral  pressure 
induced  is  proportional  to  p"f^'  •  where  p"  ■  ratio  of  the  volume  of  lateral 

reinforcement  to  the  volume  of  confined  concrete  and  f"  *  unit  stress  in 

s 

transverse  reinforcement  (which  is  assumed  to  be  equal  to  the  yield  stress)  . 

In  Figure  A. 2  are  shown  the  results  of  tests  on  rectangular  prisms 
under  concentric  load  from  reference  A. 7  and  the  modification  recommended  in 
reference  A. 5  for  members  in  flexure.  The  higher  values  in  each  of  the 
concentric  load  tests  are  for  high  strength  concretes,  the  lower  values  for 
medium  strength  concretes.  The  limits  on  the  bending  results  are  given  to 
show  the  trend  when  only  a  portion  of  the  confined  uepth  of  the  section,  d", 
is  in  compression  (c  -  depth  to  the  neutral  axis) .  A  conservative 
estimate  of  the  increase  in  strength  is  given  by  the  straight  line  whose 
equation  is: 

Af c  -  I  p"r  <  2000  psi  (A. 4) 


The  upper  limit  is  necessary  because  of  the  limited  range  of  the  test  data. 

Corresponding  to  the  increased  maximum  compressive  strength  is  a 
confined  concrete  strain,  e^',  which  can  be  expressed  as 


E0  + 


Ae 


(A. 5) 
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FIG.  A.2 :  ULTIMATE  STRENGTH  FOR  CONFINED  CONCRETE 


where  Ac  ■  increase  in  strain  at  maximum  stress  o\er  the  unconfined  value, 
c 

In  addition  to  the  volumetric  ratio,  p",  reference  A. 3  indicates  that  confined 
concrete  strain  is  dependent  on  the  ratio  of  the  minimum  dimension  of  the 
confined  core  to  the  spacing  of  the  transverse  reinforcement,  b"/s. 

Experimental  results  from  the  tests  in  reference  A. 8  and  the  recommended 
values  of  reference  A. 5  are  shown  in  Figure  A. 3  for  the  increase  in  confined 
concrete  strain.  Because  of  the  scatter  in  the  data  a  lower  hound 
straight  line  given  by  the  following  expression  is  used: 


Ac  -  0.17p"  vV'/s  <  0.008 


(A. 6) 


The  upper  limit  corresponds  to  a  total  strain  before  reaching  the  de¬ 
scending  branch  of  0.01. 

For  the  descending  branch  B'C,  the  slope  is  established  by  the  strain, 
G50C,  3t  f°r  the  con^ine<*  concrete  and  is  given  in  reference  A. A  as 


E50C 


3  +  0.002f ' 
_ c 

f?  -  1,000 

c 


+  |  P"  JbVl 


(A. 7) 


The  first  term  on  the  right  hand  side  is  identical  to  Eq.  A. 2  ,  thus 
the  second  term  represents  the  increase  in  the  50%  strain  for  the  con¬ 
fined  concrete  over  the  value  for  unconfined  concrete.  The  point  C  on  the 
descending  branch  is  determined  by  extending  a  straight  line  from  B' 
through  the  50%  point  until  the  concrete  stress  has  fallen  to  20%  of  f\ 

For  the  region  CD,  it  is  assumed  that  the  concrete  can  sustain  a  stress 
of  0.2f  'c  for  indefinitely  large  strains.  This  has  been  assumed  previously 
in  the  analysis  used  in  reference  A. 3  and  member  failure  occured  (fracture 
of  tensile  steel)  before  the  concrete  strains  became  unrealistic. 
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A. 2. 3  Cyclic  Loading  Response  of  Concrete 

The  behavior  of  concrete  under  repeated  loading  is  also  shown  in 
Figure  A.l.  Unloading  and  reloading  that  occurs  before  point  B'  (or 
point  B  in  the  case  of  unconfined  concrete)  is  assumed  to  follow  the 
Initial  tangent  slope  Ec<  Reversed  loading  on  the  descending  branch  of 
either  the  unconfined  or  confined  stress-strain  curve  is  referred  to  as 
"drop-elastic."  For  example,  on  unloading  from  point  E,  it  is  assumed 
that  0.75  of  the  previous  stress  is  lost  without  a  decrease  in  strain  (the 
"drop"  portion)  and  then  a  linear  path  of  slope  0.25  is  followed  to 
point  G  (the  "elastic"  portion).  If  the  concrete  continues  to  unload,  the 
tensile  strains  Increase  without  any  tensile  stress  developing.  On  re¬ 
loading  the  strain  must  regain  the  value  at  G  before  compressive  stress 
can  be  sustained  again.  Note  that  the  average  slope  of  the  assumed  loop 
between  E  and  G  is  parallel  to  the  initial  tangent  modulus  of  the  stress- 
strain  curve. 

This  representation  of  the  cyclic  loading  behavior  is  taken  from  reference 
A. 4.  It  can  be  modified  by  changing  the  value  of  the  slope  from  F  to  G. 

A  user  can  input  the  value  of  the  slope  as  a  constant  k  times  the  initial 
tangent  modulus.  The  value  of  k  -  0.25  is  the  default  condition. 

A.  3  STRESS-STRAIN  RELATIONSHIP  FOR  STEEL 

The  default  stress-strain  curves  utilized  for  steel  are  shown 

in  Figure  A. 4.  These  curves  cover  a  yield  point,  f^,  range  from  33  ksi  to 

75  ksi,  and  strains  from  zero  to  the  breaking  point.  They  include  two 

structural  steel  grades  with  yield  points  equal  to  33  ksi  and  36  ksi.  All 

of  the  curves  have  an  elastic  portion  AB  with  a  constant  modulus  of 

3 

elasticity,  E  -  29  x  10  ksi.  The  strain  at  the  beginning  of  yield,  e  , 
s  y 
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FIG. 


Is  equal  to  fy/Es*  The  yield  plateau  BC  varies  with  the  yield  strength  of 
the  steel.  Typical  values  for  the  strain  at  which  strain  hardening  begins, 

are  given  in  Table  1. 

The  strain  hardening  curve  CDF  reaches  a  maximum  stress,  f^,  at  a 
strain,  before  dropping  off  slightly  at  the  breaking  strain,  e . 

Typical  values  of  these  quantities  are  also  given  in  Table  1.  The 
following  expression  for  the  strain  hardening  portion  (represented  by 
a  series  of  straight  line  segments  in  the  program)  was  adapted  from  one 
developed  in  reference  A. 9 


f 

s 


(A.  8) 


Table  1 


TYPICAL  VALUES 

FOR  STEEL 

STRESS-STRAIN 

CURVES 

r,  ksi 

f  ,  ksi 
u 

ey 

esh 

cu 

cb 

33 

58 

0.0011A 

0.014 

0.15 

0.21 

36 

60 

0.00125 

0.014 

0.15 

0.20 

AO 

80 

0.00138 

0.023 

0.14 

0.20 

50 

92 

0.00173 

0.013 

0.12 

0.154 

60 

106 

0.00208 

0.0060 

0.087 

0.136 

75 

130 

0.00260 

0.0027 

0.073 

0.115 

75 


When  the  loading  is  reversed  after  the  yield  strain  has  been  reached, 
the  shape  <>f  the  stress-strain  curve  is  changed  because  it  no  longer  has  a 
well-defined  yield  point  upon  reloading.  Figure  A. 5  shows  the  general 
behavior  assumed  for  the  steel  when  reverse  loading  occurs.  On  first 
loading  to  point  1,  the  virgin  curve  described  previously  is  followed.  On 
unloading  from  point  1  to  point  2,  the  path  is  parallel  to  the  initital 
elastic  slope.  When  loading  in  the  opposite  direction  from  point  2  to  point 
3,  the  yield  point  is  missing  and  the  curve  is  described  by  Eq.  A. 8  with 
the  origin  shifted  to  point  2.  Subsequent  cycles  of  unloading  and  re¬ 
loading  follow  the  same  pattern  and  are  shown  in  Figure  A. 5. 

In  a  previous  investigation  (reference  A. 10),  a  degradation  of  stiffness 
w'ith  cycles  of  loading  was  proposed  for  the  reinforcing  steel.  However, 
a  study  of  the  original  paper  (reference  A. 11)  on  which  the  proposal  was  based 
revealed  that  tests  were  conducted  for  only  one  bar  size  (No.  11)  and 
one  yield  stress  (50  ksi) .  To  extrapolate  these  results  to  the  general 
behavior  of  all  bar  sizes  with  a  range  of  yield  points  from  33  ksi  to  75 
ksi  is  not  justified.  Furthermore,  the  data  of  reference  A.  11  showed  an 
increase  in  stiffness  in  some  cycles. 
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FIG.  A.  5s  REVERSE  LOADING  BEHAVIOR  OF  STEEL 
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NOTATION 


■  Young's  modulus  of  elasticity  for  concrete 

E  -  Young's  modulus  of  elasticity  for  steel 

8 


b" 

c 

d" 

k 

P" 

s 


breadth  of  confined  concrete  cross-section 
depth  to  neutral  axis  from  compressive  face 
depth  of  confined  concrete  cross-section 
unloading  constant  for  concrete  hysteresis  loop 
volumetric  ratio  of  transverse  reinforcement 
longitudinal  spacing  of  transverse  reinforcement 


■  compressive  stress  in  concrete 

■  compressive  strength  of  6  by  12  in.  cylinders 

-  compressive  strength  of  confined  concrete 

■  stress  in  longitudinal  reinforcement 

-  stress  in  lateral  reinforcement 

-  maximum  steel  stress  in  strain  hardening  region 
*  yield  stress  of  longitudinal  reinforcement 

=  yield  stress  of  lateral  reinforcement 

-  increase  of  concrete  strength  over  unconfined  value 


c 

.11 


'ah 


bOC 

G50u 

Ae 


■  breaking  strain  of  longitudinal  reinforcement 

-  compressive  strain  in  concrete 

-  strain  in  confined  concrete  at  maximum  stress 

■  strain  in  unconflned  concrete  at  maximum  stress 

*  strain  in  longitudinal  reinforcement 

*  steel  strain  at  onset  of  strain  hardening 

*  strain  corresponding  to  maximum  steel  stress 

-  yield  strain  of  longitudinal  reinforcement 

■  confined  concrete  strain  on  falling  branch  at  0.5f^ 

■  unconfined  concrete  strain  on  falling  branch  at  0.5f' 

c 

-  increase  of  concrete  strain  over  unconflned  value 
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APPENDIX  B 


LLEMKN I'  FAILURE  CRITERIA 

The  behavior  of  the  structural  system  Is  dependent  upon  the  behavior 
of  each  component  element  since  each  one  contributes  to  the  total  energy 
of  the  system.  The  process  of  failure  is  also  related  to  the  failure  of 
the  Individual  elements.  Some  provision  must  be  made  to  define  and  predict 
failure  In  an  element.  This  is  the  purpose  of  the  element  failure  criteria. 

B.l  DEFINITIONS 

Failure  of  a  real  concrete  member  can  be  associated  with  an  abrupt 
loss  in  its  ability  to  resist  applied  loads.  Failure  is  caused  by 
localized  behavior,  and  if  is  usually  associated  with  fracture  of  material. 

Failure  of  the  element  model  is  defined  as  those  states  which  correspond 
to  the  physical  failure  modes  of  a  real  member  and  determine  a  limit  to 
the  continuum  model  behavior.  These  states  are  detected  by  assigning  specific 
values  to  certain  quantitites  which  can  be  related  to  variables  in  the 
mathematical  model.  Since  these  variables  are  directly  associated  with 
physical  behavior,  their  values  must  be  obtained  from  test  results. 

The  derived  expressions  which  relate  these  numerical  quantities  to 
the  corresponding  quantities  determined  from  the  element  model  are 
referred  to  as  the  element  failure  criteria.  Each  failure  criterion  is 
developed  and  applied  consistent  with  the  raathemf.tical  model  of  the 
physical  system  and  the  actions. 

B . 2  CLASSIFICATION  OF  FAILURE  MODES 

A  failure  mode  is  a  distinct  manifestation  of  failure.  The  prediction 
of  each  specific  failure  mode  is  made  with  the  corresponding  failure 
criteria.  Some  criteria  can  be  directly  related  to  the  continuum  strain 
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or  stress  states  at  a  point,  while  others  rely  on  an  indirect  measure 
of  strength  through  stress  resultants  and  element  properties.  The  former 
type  of  criteria  are  referred  to  as  "micro  criteria";  the  latter  are  called 
"macro  criteria". 

The  classification  of  failure  is  made  according  to  the  dominant  stress 
state  within  an  element  at  critical  sections  defined  for  each  failure 
mode;  (critical  sections  are  discussed  in  section  B.4.2).  A  typical  element 
with  stress  resultant  -.  ’r  ation  along  the  length  is  shown  in  Fig.  B.l. 
Accordingly,  dominant  stress  states  may  be  associated  with  either  the 
bending  moment,  the  shear  force,  or  the  axial  force.  Since  normal  strain 
and  stress  values  are  defined  at  each  point  by  the  continuum  model,  the 
dominant  normal  stress  effects  (flexural  failure  and  axial  force  failure) 
are  predicted  by  micro  criteria  based  on  limiting  strain  or  stress  values. 
The  shear  stresses  are  not  predicted  in  a  direct  way;  a  nominal  (average) 
shear  stress  distribution  can  be  measured  in  an  indirect  way  based  on 
equilibrium  requirements  for  the  gross  element.  Therefore,  the  shear- 
flexure  failure  is  detected  by  a  macro  criteria. 

The  three  basic  failure  categories  are: 

1.  Flexural  failure:  dominant  normal  stress  state  caused  by 
bending,  (micro); 

2.  Shear-flexure  failure:  dominant  shear  stress  effect  in  addition 
to  the  normal  bending  stress  caused  by  a  variation  in  bending 
moment  (macro) ;  and 

3.  Axial  force  failure:  dominant  normal  stress  state  caused  by  a 
large  axial  force,  (micro). 

Since  all  three  stress  resultants  (f  ,,  f  f  .)  can  be  associated  with 

xl'  x2*  x3 

any  of  the  three  categories  of  failure,  the  distinction  between  case  3 
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FIG.  B.  1=  ELEMENT  STRESS  RESULTANTS 
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and  the  other  two  is  made  on  the  basis  of  the  strain  Btate  at  the  critical 


sections;  cases  1  and  2  are  associated  with  strain  states  which  have  a 
point  of  zero  strain  within  the  dimensions  of  the  cross  section;  case  3 
corresponds  to  the  condition  of  tension  or  compression  across  the  entire 
section,  since  the  zero  strain  point  falls  outside  the  section  dimensions. 

Each  failure  category  can  be  further  classified  according  to  the 
possible  failure  modes.  The  set  of  failure  modes  for  each  category  is 
determined  by  the  t  ype  of  reinforcement  and  the  type  of  failure  possible 
for  the  dominant  stress  conditions  prescribed.  The  sets  of  failure  modes 
for  the  failure  categories  are  defined  in  the  flow  charts  of  Figs.  B.2,3, 
and  4.  The  failure  criteria  corresponding  to  the  failure  modes  are 
developed  in  section  B.3. 

B. 3  FAILURE  CRITERIA 

Each  failure  mode  defined  in  section  B.2  requires  a  failure  criterion. 
In  addition  to  measuring  the  limiting  condition  for  the  dominant  effect, 
the  criterion  must  Include  other  effects  characteristic  of  possible  system 
behavior.  This  includes:  secondary  stress  effects,  which  provide  any 
alterations  to  the  basic  criteria  caused  by  stresses  other  than  the 
dominant  stress;  and  the  effects  of  previous  loading  history,  which 
measure  any  change  to  the  basic  form  caused  by  the  stress  variation^ 
experienced  in  an  element  during  previous  load  conditions,  such  as  load 
reversals.  In  addition,  a  user  modification  capability  is  incorporated 
so  that  the  magnitudes  of  the  basic  criteria  may  be  changed  to  allow 
certain  effects  to  be  studied.  The  secondary  stress  effects  and  the 
loading  history  effects  are  physical  measurements  obtained  from  published 
test  results.  The  user  modification  capability  is  provided  through 
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Flexural  Failure  - 


normal  stress  effect 
with  a  neutral  axis 
within  the  cross  section. 


C  -  Compression 
T  -  Tension 


Type  of 

Reinforcement  at 
Critical 
Section  / 


tens  ion''^ 
(r  compression 
failure 


tension 

or  compression 
failure 


Concrete 

Crushing 


Steel 

Fracture 


Bar  Buckling  and 
Concrete  Crushing 
Simultaneously 


Steel 

Fracture 


Fig.  B.2:  Flexural  Failure  Classification 
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Shear-Flexure  Failure  - 
combination  of  normal  and 
shear  stress  effects  with 
neutral  axis  within  the  section. 


C  **  Compression 
T  =  Tension 
DT  =  Diagonal  Tension 


^ Type  of 
Reinforcement 
at  Critical 
S\Section 


^/Tension, 
Compression,  or 
\Dlagonal  y1 
^Nlension  / 


'  Tens  1  on 
Compression,  or 
Diagonal  .Z" 
N.  Tension^ 


Cracking  Plus 
Dowel  Splitting 


Sudden  Crack 
through  Section 


Sudden  Crack 
through  Section 


Cracking  Plus 
Crushing  of 
Compression 
Zone 


Cracking  Plus 
Yielding  of 
Web 

Reinforcement 


Detectio..  of  Principal 
Diagonal  Crack 
Leading  to  Failure 


Detection  of  Principal 
Diagonal  Crack  Plus  Yielding 
of  Web  Reinforcement 
Leading  to  Failure 


Fig.  B.3:  Shear-Flexure  Failure  Classification 
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Fig.  B.4:  Axial  Force  Failure  Classification 
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coefficients  built  into  the  criteria  with  prescribed  limits  of  variation. 

All  of  the  test  data  used  in  the  development  of  the  failure  criteria 
were  obtained  from  published  sources.  Not  all  of  the  specific  information 
required  was  available;  the  deficiencies  encountered  are  noted  in  the 
criteria  developed. 

Since  the  failure  being  measured  was  usually  associated  with  some  form 
of  fracture,  the  test  data  contained  a  variable  amount  of  uncertainty 
indicated  by  some  degree  of  scatter  in  the  plotted  form.  To  represent  these 
results  in  the  form  of  a  deterministic  expression,  a  reasonable  lower 
bound  function  was  chosen  in  each  case. 

The  individual  failure  criterion  is  developed  according  to  the 
following  outline: 

1.  Basic  criterion  for  dominant  stress, 

2.  Effect  of  secondary  stresses, 

3.  Effect  of  loading  history, 

4.  User  modification  and 

5.  Assumptions. 

The  relationship  in  each  case  has  the  general  form  expressed  by: 


(Specific  criterion  value)  -  (computed  model  value.)  <  0 

The  individual  expressions  are  also  put  in  a  dimensionless  form  so  that 
units  are  not  involved  in  their  application.  The  one  exception  to  this 
rule  is  the  strain  criterion  for  concrete  crushing  in  B.3.1.1.;  the 
constants  in  this  expression  are  not  dimensionless,  even  though  the 
total  expression  is  dimensionless.  All  modification  coefficients  are 
dimensionless. 
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The  essential  details  of  all  of  the  failure  criteria  are  summarized 


in  section  B.3.4. 

B.3.1  FLEXURAL  FAILURE 

This  failure  category  is  specified  by  a  dominant  normal  stress  state 
at  a  section  in  an  element  caused  by  bending.  Another  distinction  is 
that  the  section  has  a  point  of  zero  strain  within  the  dimensions  of  the 
cross  section.  The  failure  modes  which  require  a  failure  criteria  are 
shown  in  Figure  B.2. 

B. 3.1.1  CONCRETE  CRUSHING 

This  failure  mode  is  the  crushing  of  the  concrete  in  an  unconfined 
compression  zone.  The  concrete  crushing  may  occur  progressively  from  the 
outside  surface  inward,  or  it  may  occur  in  a  sudden  disintegration  of  a 
highly  stressed  region.  Since  the  concrete  compressive  stress-strain 
response  has  a  negative  slope  beyond  an  ultimate  stress  point,  (see 
Appendix  A,  Fig.  A.I.),  a  sufficiently  large  curvature  can  cause  a 
significant  region  of  the  compression  zone  to  be  within  the  negative 
slope  influence.  At  some  point  during  an  Increasing  load,  the  bending 
moment  resistance  at  the  section  reaches  a  peak  value.  If  there  are  no 
other  regions  which  can  carry  any  additional  load  increment,  then  the 
section  will  disintegrate  suddenly.  This  condition  is  characteristic 
of  a  singly  reinforced  beam  when  a  compression  failure  occurs  before  the 
steel  yields  in  tension,  [B.34].  A  discussion  of  this  mode  of  behavior 
is  given  in  references  B.13,  18,  19  . 

Bel  ore  the  disintegration  state  is  reached,  the  outer  concrete  may 
begin  to  crush  locally.  If  this  happens,  part  of  the  moment  resisting 
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capacity  is  lost  at  that  section.  The  uniform  cross  section  property 
is  lost,  and  a  stress  concentration  effect  is  created  in  the  element.  In 
addition,  the  spalling  may  be  Irregular  causing  a  loss  in  symmetry  with 
respect  to  the  plane  of  the  structure  loading. 

If  there  are  steel  bars  in  the  compression  zone  of  the  concrete, 
crushing  can  still  occur,  but  at  a  proportionally  larger  moment  due  to  the 
load  carried  by  the  steel.  After  the  crushing  begins,  the  behavior  of 
the  steel  bars  in  compression  is  uncertain  without  web  reinforcement  to 
contain  them. 

The  entire  nonlinear  response  up  to  an  ultimate  state  can  be  predicted 
by  the  element  model  since  the  complete  stress-strain  curves  are  included. 
However,  after  crushing  of  the  outer  layers  of  concrete,  the  related 
effects  on  the  element  cannot  be  predicted  by  the  model.  Consequently, 
the  limit  to  the  continuum  model  behavior  is  associated  with  the  concrete 
crushing  state  at  a  critical  section.  The  criterion  developed  is  assumed 
to  be  valid  for  a  compression  zone  with  or  without  compression  steel. 

The  condition  of  initial  crushing  can  be  defined  by  a  maximum  strain 
value  for  unconfined  concrete.  The  stress-strain  curves  in  Figure  B.7 
show  that  the  maximum  strain  values  decrease  with  increasing  concrete 
strength  £\  This  characteristic  is  reflected  in  both  ^Qu  e20u 

strain  points  for  unconfined  concrete  defined  in  appendix  A,  Fig.  A.l, 
and  by  the  equations: 


E50u 


3.+0.002f ' 
_ c 

f'-lOOO 

c 


(B.l) 


and  e9n  *  1.8572  -  0.8572  e 

20u  50u  o 


(B.2) 
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where  ■  ultimate  concrete  cylinder  strength 

c  »  strain  at  ultimate  stress  ■  0.002 
o 

A  comparison  of  these  two  strain  points  are  shown  in  Figs.  B.5  and  B.6 

as  a  function  of  the  variable  On  the  same  figures  are  experimental 

values  of  ultimate  strain  due  to  flexure.  It  appears  that  is  a 

realistic  lower  bound  to  the  data  for  the  higher  strength  values  of 

f  (f'  >  4000. psi).  For  the  lower  values  of  f',  is  too  large: 

c  c  c  50u 

a  cut-off  value  for  the  ultimate  strain  at  0.0035  in/in  is  defined 
as  a  lower  bound  to  all  of  the  remaining  data  points. 

In  tests  of  reinforced  concrete  beams  in  reference  B.31  ,  flexural  com¬ 
pressive  strains  on  the  outer  surface  of  highly  stressed  regions  reached 
the  level  of  0.004  in/in.  before  crushing.  The  f^  values  were  in  the 
range  of  4000  -  6000  psi. 

(1)  Basic  failure  criterion:  defined  by  ultimate  strain  for  concrete 
in  compression: 

e..  -  e  <0  (B.3) 

rl  c  ■ 


where 


3  +  0.002f ' 
_ c 

f’  -  1000 


c 


=  0.0035 


e  =  concrete  compressive  strain  at  critical  section. 

The  function  is  shown  as  the  lower  limit  curve  on  Figs.  5  and  6. 

(2)  Axial  force  effect:  included  in  the  strain  state. 

(3)  Loading  history  effect:  no  measurable  difference  for  a  few  unloading¬ 

reloading  cycles,  [24,  31,  33].  Basic  stress-strain  response 
function  forms  an  envelope  to  tne  reloading  paths.  Reversal  of 
load  produces  tension  which  does  not  affect  concrete  compression 
strength. 
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FIG.  B.  5:  ULTIMATE  STRENGTH  PROPERTIES  OF  STRESS 
DISTRIBUTION  (DATA  FROM  REF.  [B.  15]). 


FIG.  B. 6:  ULTIMATE  STRAIN  FROM  TESTS  OF  REINFORCED 
CONCRETE  MEMBERS  ( DATA  FROM  REF.  [  B.  15] ). 
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CONCRETE  STRESS  (psi) 


CONCRETE  STRAIN 


FIG  B.7  FLEXURAL  STRESS- STRAIN  CURVES 
(DATA  FROM  REF  [B.  15]  ). 


0  STRAIN 


FIG.  B.8  CONCRETE  STRAIN  FUNCTIONS 
(  DATA  FROM  REF.  [B.  28]). 
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(4)  Modif icationj 


(a)  modification  by  coefficient: 


3  +  0,002  f'  „ 
— _ c  < 

f '  -  1000 

c 


CO. 0035) 

Cl.  <  C  <  1.23) 


(B.4) 


The  modified  ultimate  strain  is  defined  by  the  relationship: 
C^x(lower  limit  function).  The  upper  limit  function  is 
defined  by  ■  1.23.  (The  upper  limit  for  is  computed 
on  the  basis  of  the  evaluation  of  corresponding  to  ^Ou 
at  (maximum)  ■  8000  psi.)  The  upper  limit  function  is 
also  shown  of  Figs.  B.5  and  B.6. 

Cb)  Complete  override  is  not  possible  since  there  must  be  a 

limit  to  the  strain  in  unconfined  concrete  as  specified  in  the 
material  input  function.  The  suggested  upper  limit  for 
failure  is  defined  above. 

C5)  Assumptions: 

The  concrete  crushing  criterion  applies  to  members  with  or  without 
compression  steel. 


B.3.1.2  STEEL  FRACTURE 

It  is  possible  for  the  tensile  reinforcement  to  fracture  at  a 
critical  section  before  any  other  limit  state  is  reached  for  the  element. 
Consequently  the  fracture  of  the  longitudinal  steel  caused  by  excessive 
tensile  strain  is  considered  to  be  a  failure  mode. 

The  occurrence  of  a  bar  fracture  in  a  member  causes  significant 
changes  in  the  behavior:  it  creates  a  severe  stress  concentration  in 
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the  local  region  due  to  the  lack  of  material  continuity;  this  in  turn 
causes  additional  cracking  in  the  concrete  and  shifts  additional  stress 
to  the  remaining  steel  bars  and  the  concrete  in  compression;  if  there  are 
no  other  bars  at  the  section,  the  fracture  causes  a  complete  discontinuity 
in  the  member.  In  addition,  the  stress  distribution  in  the  fractured  bar 
varies  considerably  from  zero  at  the  break  point  to  some  tensile  value 
consistent  with  the  unbroken  bars  at  some  distance  away. 

For  a  two-dimensional  model,  all  of  the  bars  at  the  same  distance 
from  the  reference  axis  have  the  same  strain  value;  and  since  they  are 
assumed  to  have  the  same  material  properties,  the  entire  row  fractures 
at  the  limit  strain.  The  model  can  predict  the  strain  value  in  the 
nonlinear  range  up  to  the  limiting  strain  value.  After  the  fracture  point 
is  reached,  the  model  cannot  predict  the  stress  concentration  effects 
or  the  stress  distribution  in  the  fractured  bar.  Therefore,  the  limit  state 
is  defined  as  the  tensile  fracture  of  any  row  of  the  steel  reinforcing 
bars  at  the  critical  section  of  an  element. 

The  failure  condition  can  be  predicted  on  the  basis  of  a  strain  value 
at  the  point  on  the  cross  section  corresponding  to  the  bar  location.  The 
element  model  defines  this  value  directly.  Since  the  strain  is  assumed  to 
be  uniform  across  the  bar  area,  the  limit  strain  of  the  uniaxial,  stress 
strain  function  for  steel  determines  this  value.  The  numerical  values  for 
the  limit  strains  are  shown  in  appendix  A,  Fig.  A. 4. 

(1)  Basic  failure  criterion:  defined  by  a  limiting  strain  value  for  a 
steel  bar  in  tension.  (See  appendix  A,  Fig.  A. 4) 


94 


(B.5) 


where  ■  limiting  tensile  strain  value 

c  =  tensile  strain  in  longitudinal  steel  bar  at  a  critical 
8 

section, 

(2)  Axial  force  effect:  included  in  the  strain  state. 

(3)  Loading  history  effect:  it  appears  that  no  information  is  available 
on  the  question  of  the  effect  of  a  few  (2  or  3)  inelastic  stress  reversals 
on  the  fracture  strain  of  reinforcing  steel.  The  behavior  of  reinforcing 
steel  subjected  to  stress  reversal  is  documented  in  reference  B.32  .  If  the 
number  of  stress  reversals  is  small,  simplifications  may  be  introduced  into 
the  stress-strain  function  which  includes  stress-reversal,  basing  the 
stress  reversal  curve  on  the  original  monotonic  stress-strain  response. 

Others  have  utilized  this  concept  to  define  an  idealized  response  for 

steel  with  stress  reversal,  [B.  1,8, 29].  But  none  of  these  have  indicated 
the  fracture  strength.  For  simplification,  it  is  assumed  that  the  strain 
at  fracture  remains  the  same  ar  the  monotonic  fracture  point,  regardless 
of  the  history  of  loading. 

(4)  Modification: 

(a)  no  parameter  modification  is  necessary  since  the  fracture 
point  defined  by  the  input  function  for  the  material  is  net 
altered  in  the  failure  criteria 

(b)  no  override  is  possible  because  the  physical  limit  of  the 
stress-strain  function  for  each  material  is  defined  independent 
of  the  limit  conditions  for  the  model. 

(5)  Assumptions: 

The  limiting  tensile  strain  value  for  the  monotonic  stress-strain 
function  for  steel  is  a  valid  measure  of  the  fracture  strength  including 
loading  history  effects. 


Note: 


No  special  check  is  made  for  the  case  where  bending  in  a  singly 
reinforced  concrete  element  occurs  opposite  to  its  reinforced  strength. 
Clearly  the  concrete  will  crack  in  tension  and  form  a  discontinuity  in 
the  members  at  the  section.  Physically,  the  longitudinal  bars  could  be 
positioned  at  any  depth  in  a  member.  Whether  or  not  there  is  sufficient 
resisting  strength  in  the  couple  that  is  formed  by  the  resulting  concrete 
compression  zone  and  the  steel  in  tension  depends  on  the  moment  to  be 
resisted  and  the  steel  location  at  a  section.  If  the  steel  is  near  the 
tension  surface,  a  larger  moment  can  be  resisted  than  for  the  case  where 
the  steel  is  near  the  compression  surface.  In  any  case,  the  limitations 
provided  by  the  concrete  crushing  (B.3.1.1)  and  steel  fracture  (B.3.1.2) 
are  sufficient  failure  checks  for  any  placement  of  the  longitudinal  steel. 

B.3.1.3  BAR  BUCKLING  AND  CONCRETE  CRUSHING  SIMULTANEOUSLY 

This  failure  mode  is  the  process  of  the  bending  out  of  reinforcing 
bars  and  the  simultaneous  crushing  of  concrete  causing  a  sudden  destruction 
of  the  compression  zone.  It  is  characteristic  of  a  member  with  web  rein¬ 
forcement  (see  Fig.  B.9). 

This  behavior  has  several  distinguishing  traits: 

a.  The  condition  is  not  possible  until  the  concrete  outside  of  the 
compression  bars  has  begun  to  crush;  otherwise  the  bars  are 
adequately  restrained  from  buckling; 

b.  The  process  is  most  likely  initiated  by  the  expansion  of  the 
confined  concrete  near  the  ultimate  stress  level,  effectively 
pushing  the  bars  outward  from  their  normal  positions; 
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c.  The  bars  are  in  a  yielded  ..tate  in  compression,  usually  in  the 
strain  hardening  range,  at  the  time  of  buckling,  [B.31]; 

d.  Normal  web  reinforcement  provides  a  restraint  to  displacement 

of  the  longitudinal  bars  at  their  contact  points,  thus  providing 
a  significant  influence  on  the  bending  strength  of  a  bar  segment, 
[B.5,31] . 

The  concrete  expansion  effect  can  be  measured  by  its  volumetric  strain 
under  compressive  stress.  The  lateral  strain  and  volumetric  strain  of 
unconfined  uniaxially  loaded  concrete  begin  altered  behavior  at  approximately 
80%  of  ultimate  strength  to  the  extent  that  near  the  ultimate  stress  the 
volumetric  strain  has  changed  signs  from  compression  to  tension,  (see 
Fig.  B.8).  This  indicates  that  there  is  an  expansion  effect  developed 
to  push  the  bars  out  of  line,  [B.28]. 

An  important  effect  is  the  combination  of  yielded  steel  and  crushed 
concrete  cover  which  allows  the  failure  process  to  occur.  To  see  the 
possibility  of  this  combined  effect,  the  steel  strains  at  strain  hardening 
are  compared  to  the  range  of  strain  at  which  unconfined  concrete  crushing 
occurs : 

Steel  Yield  Strength 

75  60  50  40 

Steel  strain  k  .  k  ,  k  .  k  . 

si  si  si  si 

at  strain  hardening:  ^  i _ i  _ | _  i _ ^ 

Strain: 

(in/in)  0.00  0.01  0.02  0.03 

- - 1 —  - 1 - J - 1 — 

Concrete  strain 

at  crushing  , 

(unconfined)  2000.^^8000. 

If  some  allowance  is  made  for  the  fact  that  the  normal  strain  at  a  section 
is  larger  at  the  surface  than  at  the  bar  location,  it  is  reasonable  to 
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state  that  when  the  steel  enters  the  strain  hardening  range,  the  con¬ 
crete  cover  has  reached  the  crushing  state. 

In  reference  B.5,  experiments  were  performed  to  study  the  behavior  of 
the  longitudinal  steel  as  affected  by  compressive  stress  and  lateral 
reinforcement.  One  conclusion  reached  was  that  the  size  of  the  lateral 
steel  bar  was  not  important  for  providing  restraint  against  the  outward 
displacement  of  the  longitudinal  bars;  a  positive  direct  tensile  connection 
as  shown  in  Fig.  B.10  for  a  corner  bar  and  an  interior  bar  provide  adequate 
restraint  with  two  smallest  diameter  lateral  bars. 

In  observations  of  this  failure  form  in  tests,  it  is  difficult  to 
determine  whether  the  concrete  crushing  or  the  steel  buckling  initiates 
the  final  destruction  [ B . 31 ] . 

The  element  model  can  predict  the  flexural  behavior  of  an  element 
into  the  nonlinear  material  range.  However,  it  cannot  predict  the  behavior 
after  the  buckling  of  the  compression  bars  at  a  critical  section  since 
the  compression  zone  has  been  destroyed  in  the  region  of  buckling.  Not 
only  has  a  partial  discontinuity  been  created  at  the  section  but  the 
material  uniformity  along  the  length  has  been  altered.  The  limit  state 
is  defined  to  be  the  buckling  of  compression  bars  stressed  into  the 
strain  hardening  range  at  the  critical  section  of  an  element.  The  buckling 
condition  for  the  outer  layer  of  bars  is  considered  sufficient  for  defining 
the  limit  state  if  there  are  multiple  layers  involved  since  an  unknown 
stress  concentration  and  redistribution  is  created  due  to  this  localized 
effect.  For  the  member  with  web  reinforcement,  the  crushing  of  the  outer 
concrete  is  not  considered  to  alter  the  behavior  significantly.  The 
element  model  is  assumed  to  be  valiu  up  to  the  point  of  bar  buckling. 
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The  expression  for  predicting  the  critical  stress  for  a  uniform  bar 
in  simple  compression  is  used  to  measure  the  buckling  condition.  To 
include  the  effect  of  the  strain  hardening  state  of  the  material,  the 
tangent  modulus  property  is  Included.  The  critical  stress  is  defined  by 
the  following  expression: 


cr 


C„  •  ir 


Et 

CS/K)2 


(B.6) 


where 


f  -  critical  compressive  stress  in  longitudinal  reinforcing 
bar 


-  tangent  modulus  for  steel  at  fcr 
S  -  spacing  of  web  reinforcement 
K  -  radius  of  gyration  of  the  bar 
0>2  “  end  restraint  coefficient . 

■  1:  pinned  end  condition 
C2  ■  A:  fixed  end  condition 

This  equation  for  f  was  used  in  reference  B.5  as  a  basis  for  determining  the 
spacing  of  web  reinforcement  in  compression  members;  a  value  of  C2  *  2  was 
used. 


To  Implement  this  prediction,  the  model  can  compute  the  stress  in  a 
bar  directly  from  a  given  strain  state,  and  this  value  can  be  compared  to 
a  specific  critical  stress  value.  One  bar  size  and  one  material  are 
representative  of  the  entire  layer  checked. 

(1)  Basic  failure  criterion:  the  relationship  is  based  on  the  equation 
for  critical  stress  of  an  initially  straight  uniform  bar  in  simple  com¬ 
pression  : 
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f  <  0 
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where 

2  Et 

i  *■  C„tt  *  -----  defined  in  equation  B.6; 

cr  2  (S/K) 1 

f  *  compressive  stress  in  longitudinal  reinforcing  bar. 

8 

If  f  is  expressed  in  terms  of  the  bar  diameter  D,  i.e. 
cr 

let  k  “  D/4,  and  if  the  criterion  is  put  in  a  dimensionless  form,  the 
expression  can  be  re-written  as: 


- 

16  \S) 


-  1  <  0 


(B.8) 


(2)  Axial  force  effect:  included  in  the  strain  state. 

(3)  Loading  history  effect:  since  the  stress-strain  response  is 
uniquely  defined  for  unloading  and  reloading  of  steel  bars,  the 
check  can  be  made  at  any  point  where  a  tangent  modulus  exists. 

The  check  is  applied  for  a  compressive  stress  state  in  the  strain 
hardening  range  for  any  cycle  of  loading.  This  is  demonstrated 
in  Fig.  B.ll. 

(4)  Modification: 

a.  Modification  by  coefficient  is  made  through  the  constant  C£. 
This  reflects  the  effects  of  various  end  restraints  on  the 
critical  buckling  stress.  The  end  restraint  on  the  segment 
between  web  bars  is  a  function  of  the  continuity  of  the  bar 
and  its  freedom  to  deform.  The  modified  equation  is: 

2  2  E 

c2  •  is  {I)  •  r  - 1 : 0  <B-9> 

s 

a.  <  c2  <  4.) 
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FIG.  B.  11=  FAILURE  CHECK  FOR  BAR  BUCKLING. 
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C9  -  2  is  used  as  the  lower  bound  value. 

b.  The  user  may  choose  to  completely  suppress  this  particular  failure 
mode.  By  overriding  this  condition  the  member  will  maintain  the 
original  configuration  of  the  steel  for  all  stress-strain  values. 

5.  Assumptions: 

a.  The  steel  bar  is  initially  straight  and  uniformly  compressed; 

fc.  The  steel  bar  cannot  buckle  until  it  is  stressed  into  the  strain 
hardening  range  of  behnvior; 

c.  Ali  normal  web  reinforcement  bars  provide  sufficient  restraint 
to  prevent  lateral  displacement  of  the  longitudinal  bars  at  the 
point  of  contact; 

d.  Buckling  of  the  longitudinal  steel  bars  initiates  the  destruction 
of  the  comnression  zone  at  the  critical  section. 

B.3.1.4  STEEL  FRACTURE 

This  failure  mode  is  the  fracture  of  the  longitudinal  steel  caused  by 
excessive  tensile  strain.  Since  this  fracture  is  unaffected  by  web  rein¬ 
forcement  in  a  member,  the  failure  criterion  is  identical  to  the  criterion 
developed  in  section  B.3.1.2. 

B.3.2  SHEAR- FLEXURE  FAILURE 

This  category  of  failure  defines  the  failure  modes  and  criteria  fo." 
the  effects  of  shear  stress  and  normal  stress  acting  in  combination  within 
a  member.  The  failure  modes  which  require  a  failure  criteria  are  shown 
in  Fig.  B.3. 

The  distinct  form  of  behavior  for  this  type  of  failure  is  a  diagonal 
crack  which  more  or  less  follows  a  path  normal  to  the  principal  tensile 


stress  trajectories  in  the  member,  Although  several  similar  cracks  may  be 
Initiated,  one  crack  becomes  dominant  as  the  load  is  Increased.  If  web 
reinforcement  exists  at  the  crack  location,  sudden  material  separation  is 
prevented. 

The  principal  diagonal  crack  altera  the  member  behavior  by  Introducing 
a  discontinuity.  The  crack  is  deep  enough  to  disrupt  the  flexural 
characteristics  within  the  element  by  reducing  the  size  of  the  compression 
zone  in  a  localized  region.  As  a  result,  a  redistribution  of  stresses  is 
required  to  maintain  equilibrium  in  the  altered  state,  [B. 10, 21, 35] . 

The  discontinuous  behavior  of  diagonal  tension  cracking  in  reinforced 
concrete  members  is  associated  with  the  weakness  of  concrete  to  principal 
tensile  stress.  The  combination  of  normal  tensile  stress  due  to  flexure, 
and  shear  stress  due  to  variation  in  flexure  can  produce  a  more  severe 
principal  tensile  state  at  a  point  than  for  conditions  due  to  flexure 
alone.  The  tensile  strength  of  concrete  depends  on  local  properties  near 
the  point  of  fracture  in  addition  to  the  strain  distribution,  hence  it 
is  random  in  nature.  This  manifests  itself  in  the  unpredictability  of 
specific  locations  and  shapes  of  diagonal  cracks. 

This  is  in  contrast  to  the  more  predictable  behavior  in  bending  of 
the  same  type  of  member.  In  pure  bending,  the  member  strength  depends  on 
the  compression  zone  properties  of  concrete  and  the  tensile  strength  of 
the  longitudinal  steel.  The  resulting  crack  formation  in  the  tensile  zone 
is  of  only  secondary  significance  because  these  cracks  are  stabilized  by 
a  uniform  compression  zone. 


7 


1 


B. 3.2.1,  DETECTION  OF  PRINCIPAL  DIAGONAL  CRACK  LEADING  TO  FAILURE 

If  a  diagonal  crack  occurs  in  a  member  without  web  reinforcement, 
there  are  three  forms  of  behavior  which  describe  the  ultimate  failure 
condition : 

a.  The  crack  may  propagate  through  the  entire  section; 

b.  The  crack  may  stop  near  a  three  dimensional  compression  zone,  e.g., 
near  a  concentrated  force  or  a  support  point;  the  final  failure 
form  is  a  destruction  of  the  reduced  compression  zone  above  the 
crack; 

c.  The  crack  may  propagate  parallel  to  the  tensile  steel  bars 
effectively  removing  their  contribution  from  the  element  load 
resistance. 

The  ultimate  failure  conditions  described  by  the  crushing  of  the  com¬ 
pression  zone  above  the  crack  (b)  and  splitting  along  the  longitudinal 
steel  (c)  are  considered  post-cracking  behavior  since  they  occur  after  the 
formation  of  the  principal  crack. 

The  significant  factors  which  affect  the  shear  cracking  behavior 
have  been  determined  principally  from  test  results.  Those  factors  con¬ 
sidered  most  significant  are: 

a.  the  longitudinal  steel  percentage, 

b.  the  tensile  strength  of  concrete, 

c.  the  dimensions  of  the  cross  section  in  comparison  to  the  length 
of  the  member, 

d.  an  axial  force  effect  in  addition  to  effects  of  lateral  forces 
applied, 

e*  the  location  of  the  region  of  maximum  shear  force  with  bending, 


i. 
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f.  the  anchorage  and  bond  characteristics  of  the  longitudinal  steel, 
and 

g.  the  wav  loads  are  applied  to  a  member,  directly  or  indirectly. 

These  factors  can  be  divided  into  three  main  groups:  the  initial  shape 
and  reinforcement  of  the  member  (a.,  c.),  material  properties  (b.,  f.), 
and  the  internal  stress  distribution  (d.,  e.,  g.).  The  steel  properties 
are  not  significant  because  the  steel  does  not  yield  before  the 
formation  of  the  significant  crack  unless  the  flexural  effects  are  dominant 
[B.101.  The  anchorage  and  bond  characteristics  of  the  longitudinal  steel 
have  the  most  Lnfluerce  on  the  post-cracking  behavior. 

The  limitations  Imposed  on  a  continuum  model  are  based  on  physical 
conditions  which  significantly  alter  the  dominant  normal  stress  behavior 
predicted  by  the  model.  The  discontinuous  effects  of  a  diagonal  crack  and 
the  resulting  stress  concentration  at  the  tip  of  the  crack  are  not  included. 
Therefore,  the  limit  state  is  considered  to  be  the  initial  formation  of  the 
principal  diagonal  crack.  This  precludes  the  consideration  of  the  post¬ 
cracking  behavior  as  limit  states. 

The  element  model  measures  the  normal  stress  distribution  along  the 
length  of  the  element.  Flexural  and  axial  stress  effects  are  both  derived 
from  the  given  strain  state.  To  get  a  measure  of  the  shear  stress  effects, 
only  indirect  procedures  are  available:  i.e.,the  shear  force  can  be  com¬ 
puted  to  satisfy  equilibrium  for  the  normal  stress  resultants,  and  a 
nominal  measure  of  the  stresses  can  be  based  on  an  average  distribution 
equivalent  ♦•a  a  c hear  force  at  a  cross  section.  Consequently,  the  con¬ 
ditions  causing  the  formation  of  a  diagonal  crack  cannot  depend  upon 
stresses  at  a  point  (micro  criterion) ;  the  true  stress  distribution  is 
unknown.  The  crack  prediction  must  be  based  on  the  nominal  conditions 
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indirectly  related  to  the  stresses  predicted  by  the  model. 

The  indirect  measure  usually  appears  as  an  empirical  relationship 
between  gross  element  properties  and  nominal  stress  values.  From  an 
ui-ier standing  of  the  uncertainties  associated  with  the  formation  of  the 
principal  diagonal  crack,  it  is  understandable  that  the  relationships 
developed  for  their  prediction  are  inaccurate  for  generalized  conditions. 

The  criterion  developed  in  this  section  incorporates  the  variables  observed 
to  be  significant  in  tests,  while  accounting  for  the  uncertainties  by 
providing  lower  bound  functions.. 

There  are  two  important  considerations  in  defining  the  failure  criterion 
for  diagonal  tension  cracking:  (1)  the  equation  for  predicting  the  crackin;; 
strength  of  a  member  at  a  particular  section  for  generalized  conditions;  and 
(2)  the  choice  of  a  critical  section  at  which  to  apply  the  equation.  The 
prediction  equation  is  discussed  below;  the  application  to  a  critical  section 
is  discussed  in  the  implementation  section  (B.4). 

The  form  for  the  failure  criterion  is 


V  (x)  v  (x) 
c  _  c 

bd/F  SF 

c  c 


A  +  B 

/V 


V(x) 

M(x) 


(B.10) 


v 

c 

V 

c 

b 

A,B 

f' 

c 


nominal  measure  of  cracking  stress  at  station  x 

cracking  force  at  a  section  x 

gross  cross  section  width 

constants  to  be  defined  by  data 

ultimate  cylinder  strength  of  concrete 

is  a  measure  of  tensile  strength) 

longtudinal  steel  percentage  ■  A  /bd 

6 


p 
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V(x),  M(x)  -  shear  force  and  bending  moment  at  a  section  x 
where  is  computed 

d  -  effective  depth  of  longitudinal  steel. 

The  form  of  this  equation  is  identical  to  that  presented  in  reference  B.30  , 
and  it  has  also  been  used  in  modified  versions  by  others  [B.10,  17]. 

(Note:  other  equations  have  been  developed  for  computing  v^,  e.g., 

references  B.21,37  ,  but  the  fit  of  data  is  no  better  in  general.)  This 
equation  defines  the  nominal  shear  stress  at  the  initial  crack  formation 
for  a  specific  location  x  as  a  function  of  the  element  dimensions , 
reinforcement,  and  material  properties  (p,  b,  d,  f ’) ,  and  as  a  function  of 
the  Internal  stress  distribution  (V(x)/M(x))  evaluated  at  the  same  location 
x.  The  ratio  (V/M)  measures  the  effect  of  the  shear  force  and  the  bending 
moment  on  the  principal  tensile  stress,  and  hence  on  the  crack  formation. 

As  a  ratio,  this  Influence  is  translated  into  a  distance  measurement, 
from  the  point  where  M  *  0  to  the  crack  location. 

In  order  to  evaluate  the  constants  A  and  B  consistent  with  these 
characteristics,  the  V(x)  and  M(x)  values  must  be  known  at  the  crack 
location,  in  addition  to  the  other  factors.  The  two  basic  parameters, 

v  W  a  \ 

F.(x)  -  -  and  F2(x)  -  1000  • 

1  bd  /F  2  /f^  M(x) 

c  c 

can  be  evaluated  for  specific  test  cases  for  a  wide  variation  in  each 
factor.  The  choice  of  a  specific  function  to  fit  the  data  will  establish 
A  and  B. 

The  process  chosen  here  is  to  develop  a  basic  equation  using  only  the 
results  from  simple  beam  tests  with  concentrated  loads  since  they  represent, 
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in  some  form,  all  concentrated  load  conditions.  In  addition,  there  are 
many  test  results  of  this  form  available.  This  basic  equation  can  then 
be  checked  against  other  conditions  for  more  general  applicability.  The 
basic  data  is  from  reference  B.30  and  is  plotted  in  Fig.  R.12.  The  lower 
bound  function  is  shown  on  the  figure  and  can  be  expressed  as: 

F1(x)  -  1.5  +  3.5  F2(x)  <  3.0  (B.ll) 

where  A  =  1.5  and  B  =  3.5  define  the  intercept  distance  and  the  slope  of 

the  function. 

The  method  of  load  application  is  significant  to  the  cracking 
strength  particularly  for  small  (a/d)  values,  where  a  is  the  shear  span 

length,  [B.12],  Such  refinement  is  not  considered  appropriate  in  this 

behavior  prediction  for  lack  of  quantitative  measurements  and  precise 
characterization  of  the  physical  loads  actually  applied.  It  should  be 
noted,  however,  that  for  most  beam  dimensions  and  loading  O^OO^.S), 
the  (a/d)  values  are  larger,  and  the  support  and  load  effects  are  not 
significant  in  the  region  of  the  crack. 

The  check  of  the  lower  bound  function  for  the  cases  of  simple  beams 
wlLh  uniformly  distributed  loads  and  for  restrained  and  continuous  beams 
with  concentrated  loads  are  shown  separately  in  Figs.  B.13  and  B.14, 
with  reasonable  acceptance. 

(1)  Basic  failure  criterion:  defined  by  the  lower  bound  function 
developed: 

F^x)  -  I (x)  <  0  (B.12) 
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FIG.  B.  12=  DATA  FOR  INITIAL  SHEAR  CRACKING -(A) 
(PLOTTED  DATA  FROM  REF.  [B.  30]  TABLES 
5-1,2,8,9,12.) 
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FIG.  B.  I4>  DATA  FOR  INITIAL  SHEAR  CRACKING-(C) 
(PLOTTED  DATA  FROM  REF.  [B. 30]  TABLES 
5-3,4.) 
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where 


F^x)  -  1.5  +  3.5  •  F2(x)  <  3.0 
F(x)  =*  V/btl/F 

C  • 

F(x)  is  a  nominal  measure  of  the  actual  shear  force  at  the  critical 
section. 

The  equation  is  considered  applicable  to  any  concentrated  load 
condition  for  straight  elements  because  the  shear  and  moment  functions 
can  be  broken  down  into  equivalent  simple  beam  segments  when  the  sections 
of  zero  moment  are  located.  Each  segment  of  constant  shear  can  be 
examined  separately. 

(2)  Axial  for.e  effect:  this  effect  is  to  measure  the  change  in  the 
basic  criterion  caused  by  the  application  of  a  tensile  or  compressive 
axial  force.  The  effect  is  reflected  in  the  alteration  of  the  principal 
stress  trajectories  in  the  element,  A  compressive  force  causes  a  larger 
compression  zone  which  forces  the  diagonal  crack  to  bend  over  at  a 
shallower  depth  from  the  tension  face.  A  tensile  force  causes  the 
compression  zone  to  be  smaller  in  depth  resulting  in  a  deeper  and  less 
inclined  crack  [B.10,  .14] . 

The  axial  force  influence  is  not  symmetric,  particularly  for  larger 
values.  In  the  limiting  case,  where  the  axial  force  dominates  the 
behavior,  the  effects  of  compression  and  tension  are  distinctly  different. 
(This  behavior  is  discussed  in  section  B.3.3.). 

The  choice  for  an  equation  to  predict  the  axial  force  affect  was 
made  on  the  basis  of  the  following  points: 

a.  It  is  easier  and  more  efficient  to  work  with  strains  and  stresses 


from  the  model  rather  than  stress  resultants. 


b.  The  equation  should  be  numerically  easy  to  apply  and  easy  to 
modify  by  coefficients. 

c.  Prediction  of  the  axial  force  effects  should  be  no  less  accurate 
in  general  application  than  the  currently  accepted  design 
equations. 

The  two  alternative  approaches  to  specifying  a  prediction  equation 
were:  (1)  to  use  the  current  ACI  Code  equation  (ACI  Standard  318-71), 


V  -  1.9/f7  +  2500  p.  —  -tt-.  .. 


(B.13) 


or  (2)  to  develop  a  different  equation  and  check  it  with  code  equation. 
The  choice  was  made  to  develop  a  different  equation  for  the  following 


reasons : 


a.  The  ACI  Code  equation  required  the  computation  of  stress 
resultants  with  no  alternative; 

b.  To  modify  the  equation  with  a  coefficient  for  the  effect  of  N, 
and  to  determine  lower  and  upper  limits  to  its  variation  is 
difficult  in  the  form  presented  by  the  Code  equation. 

c.  A  different  equation  could  be  developed  to  more  closely  fit 
the  desired  qualities; 

d.  There  appeared  to  be  sufficient  published  test  results  as  a 
basis  for  the  equation  development  and  for  testing  its  accuracy 
for  both  tension  and  compression  effects. 

The  initial  form  of  the  equation  was  based  on  a  form  suggested  in 


reference  3.10  ,  i.e., 


Vcn-(X)  -  [1.5  +  3.5F,(x)]-[l-a^X-}-J  (B.14) 

bd/F  2  VcnU) 

c 
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The  new  terms  ate  defined  to  be: 

V  (x)  *  cracking  shear  force  at  section  x  with  axial  force 
cn 

N(x)  =  axial  force  (tension  position) 

a  =  coefficient  for  axial  force  effect. 

This  equation  was  modified  because  N/Vcn  did  not  correctly  measure  a 
change  in  the  cracking  strength  for  all  cases.  This  is  due  to  the  fact 
that  Vcn  is  influenced  by  N;  i.e.,Vcn  increases  for  N  compression, 
and  decreases  for  N  in  tension.  This  equation  is  more  reasonable  if 
there  is  a  fixed  ratio  for  N/V  up  to  the  cracking  strength. 

The  modified  form  uses  V£(x)  as  the  reference  force  for  N(x), 

where  V  (x)  is  the  cracking  force  without  an  axial  force  effect;  i.e. , 

c 


V 

cn 


(x) 


bd/F 

c 


[1.5  +  3.5  F2(x)]  *  [1. 


N(x)  , 

n  V  (x)  ] 
c 


(B.15) 


Noting  that: 

V  (x) 

F  (x)  = - 3  =  1.5  +  3.5  F9(x)  (equation  B.ll) 

bd/f ' 

c 


and  nondimensionalizing  N  and  by  the  factor  bd/F,  equation  (B.15)  can 
be  written: 


V  (x)  V  (x)  N(x)/bd/fr 

-  =  _c -  •  [i.  -  a  •  - c  j 

bd/F  bd/F  V  (x)/bd/F 

c  C  c  c 


and  finally, 


v  (x) 

cn 

bd/F 


c 


V  (x) 

C 

bd/F 


.  N(x) 

bd/F 


(B.16) 
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and  F^(x) 


V  (x) 

For  F.  (x)  -  — - 

3  bd/F 

c 


N(x) 

bd/F" 


equation  (B.16)  can  be  written: 

F3(x)  =  F1(x)  -  a  F4(x) 


(B. 17) 


The  limitation  (£  3.0)  shown  in  equation  B.ll  is  omitted  in  the  derivation, 
but  it  is  to  be  incorporated  in  the  application. 

To  evaluate  the  lower  limit  for  the  coefficient  a,  the  equation  B.17 
was  rearranged  as  shown  below: 

v  -  v 
cn  c 

-  *  a 

vc 

V  V 

cn  c 

where  v  -  — —  and  v  =  rj  .  These  two  basic  factors  are  evaluated  for 

cn  bd  c  bd 

a  series  of  tests  and  plotted  in  Fig.  B.15.  The  lower  limit  function  is 
shown  as  a  *  0.10  for  both  tension  and  compression.  A  reasonable  average 
value  was  chosen  for  the  lower  limit  of  a  to  avoid  compounding  lower 
bound  features  into  one  equation.  The  basic  equation  for  Vc  was  previously 
defined  to  be  a  lower  bound  fit  of  data. 

The  f*nal  justification  for  the  implementation  of  the  derived 
equation  is  based  on  its  ability  to  predict  the  observed  trends  in  test 
results.  The  derived  equation  was  applied  to  the  data  in  reference  B.14 
which  included  tests  for  tension,  compression  and  zero  axial  force.  The 
results  are  shown  in  Fig.  B.20  as  a  comparison  with  test  values;  on  the 
same  plot,  the  ACI  Code  equation  is  applied  to  the  same  set  of  data.  It 
is  clear  that  the  derived  equation  is  as  valid  as  the  Code  equation  for 
measuring  all  three  axial  force  effects.  The  data  from  reference  B.14  is 
shown  to  be  incorporated  into  the  derivation  of  the  coefficient  a  and 


N_ 

V 


(B.18) 
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FIG.  B.  15=  DATA  FOR  INITIAL  SHEAR  CRACKING  WITH  AXIAL 

FORCE  EFFECTS.  (DATA  FROM  REF.  [B.  3,  14, 10,  23,  27] ) 
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in  checking  the  equation.  However,  the  data  points  on  Fig.  B.15 
corresponding  to  reference  R.14  were  plotted  after  the  equation  check  was 
made  so  that  an  independent  set  of  tests  were  ivailable  to  make  the 
check. 

In  reference  B.14  ,  a  recommendation  is  nade  to  incorporate  a  cut-off 

value  for  large  tensile  forces;  i.e.  for  N  (tension)  >  4 fV  ,  V  is  to 

c  c 

be  taken  equal  to  zero.  It  is  argued  that  too  much  preliminary  cracking 
caused  by  the  tensile  force  destroys  the  shear  resistance  of  the  member. 
There  is  no  such  cut-off  recommended  with  the  failure  criteria  for  the 
following  reasons: 

a.  All  of  the  test  results  from  reference  B.14  do  not  support  the 
recommeudatii'  made; 

b.  The  recommendation  is  based  on  one  test  for  which  the  axial 
load  was  applied  first,  then  the  lateral  load  up  to  cracking, 
which  is  not  realistic  in  an  actual  s;ructure; 

c.  The  plotted  results  of  the  derived  eq lation  in  Figure.  B.20 
did  not  indicate  any  trend  to  support  the  cut-off  limit. 

The  basic  failure  criteria  with  axial  force  effect  is: 

F3(x)  -  F(x)  <  0  (B. 19) 

where  Fj(x)  *  F  (x)  -  «  •  F^(x) 

F(x)  =  V(x)/bd/fT. 

c 

(3)  Loading  history  effects;  the  effect  of  reversal  of  load  on  the  shear 
strength  of  members  was  tested  and  results  presented  in  reference  B.2  .  ""wo 
conclusions  are  significant; 
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It  was  found  that  cracking  load  associated  with  a  reverse  load  is 


of  the  same  order  of  magnitude  as  the  initial  cracking  load,  even 
in  the  presence  of  the  initial  diagonal  crack; 

b.  The  repetition  of  a  few  high  reverse  loadings  does  not  cause  a 
significant  decrease  in  the  strength  of  beams  failing  in  shear 
as  compared  to  monotonic  loading  failure. 

Therefore,  it  is  assumed  that  the  developed  criteria  is  valid  for  any 
loading  history  consistent  with  a  few  cycles  of  loading  effects. 

(4)  Modification: 

00  Modification  by  coefficient:  the  suggested  modification  is  made 
by  increasing  the  slope  of  the  basic  function,  equation  B.ll, 
but  not  the  point  of  intersection  at  F^(x)  =  1.5.  The  horizontal 
cut-off  line  can  be  defined  by  the  intersection  with  F2(x)=3/7 
for  each  slope.  The  form  for  modification  is  defined  below: 

F*(x)  =  ([1.5  +  C3  ’  F2(x)]  <  1{1^}  j  -  o-F4(x)  (B.20) 

and  a  (tension)  =  =  0.025(4.-C^) 

a  (compression)  =  ra  =  0. 050(2. +C^) 

The  range  for  coefficients  and  are  defined  by: 

(3.5  <  C3  <  ^  ) 

(0.  <  C4  <  1.) 

The  upper  limit  function  for  N(x)  =  0  is  shown  on  Fig.  B.12 
corresponding  to  =  35. /6.;  the  upper  limit  for  and  are 
shown  on  Fig.  B.15  for  =  1.0. 

Cb)  This  failure  criterion  can  be  overridden  since  it  is  an  indirect 
measure  of  a  possible  discontinuity.  If  it  is  ignored,  the  user 
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should  realize  that  only  the  limitations  of  the  flexural 
behavior  are  in  effect,  and  unrealistic  behavior  is  possible. 

(5)  Assumptions: 

The  shear  crack  in  a  member  with  several  shear  spans  (see  Fig.  B.18), 
can  be  predicted  by  a  criteria  based  on  simple  beam  tests  with  or  without  an 
axial  force. 
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B.3.2.2  DETECTION  OF  PRINCIPAL  DIAGONAL  CRACK  PLUS  YTELDTNC  Ob'  WEB 

ri:inforci:mi:nt  ij-adinc  to  faiijiki:. 

IF  ;i  diagonal  crack  occurs  in  a  member  with  web  reinforcement,  material 
separation  is  prevented  bv  the  bars  intersected  by  the  crack.  The  only  form 
of  web  reinforcement  considered  is  the  closed  form  normal  to  the  longitudinal 
reinforcement . 

There  are  two  behavior  conditions  which  describe  ultimate  failure  states 
for  a  member  with  web  reinforcement  f B . 14 ,  30]: 

a.  The  diagonal  crack  may  propagate  through  the  member  and  cause  sudden 

yielding  of  the  web  reinforcement  at  the  crack  location; 

b.  The  diagonal  crack  may  be  contained  by  the  web  reinforcement  until 
crushing  of  the  reduced  compression  zone  at  an  increased  load 
before  or  after  the  web  steel  yields. 

In  addition,  experimental  tests  have  shown  that  the  stress  in  the  web  bars 
at  the  crack  remain  very  small  until  the  significant  diagonal  crack  forms 
[B.7,  9,  11,  25],  This  is  due  to  the  fact  that  brittle  materials  possess 
low  extensibility  ( B . 3 6 ] . 

To  predict  shear-flexure  failure  of  a  member  with  web  reinforcement 
requires  some  measure  of  the  behavior  after  the  formation  of  the  diagonal 
crack.  This  means  that  the  stress  in  the  web  bars  intersected  by  the  crack 
must  be  predicted  since  subsequent  behavior  depends  on  the  containment  of 
the  crack  by  these  bars.  However,  to  predict  the  destruction  of  the  com¬ 
pression  zone  above  the  crack  requires  some  knowledge  of  the  size  of  the 
compression  zone  and  the  stress  concentration  effects  in  the  localized 
region.  This  is  beyond  the  model  capability.  Therefore,  the  limit  state 
considered  for  this  failure  mode  is  the  detection  of  the  yielding  of  web 
reinforcement  in  the  range  of  the  crack.  As  long  as  the  web  bars  are  below 
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the  yield  stress,  the  element  is  considered  to  be  reliable.  Beyond  the 
yielded  state,  the  element  model  is  considered  to  be  invalid,  and  the 
behavior  cannot  be  accurately  predicted. 


An  approximate  measure  of  the  web  steel  stress  can  be  computed 
indirectly.  The  form  of  a  prediction  expression  can  be  developed  on  the 
basis  of  measured  stirrup  strain  during  loading  (See  Fig.  B.16). 
Accordingly,  the  assumed  stress  variation  is  zero  up  to  the  cracking  load, 
and  thereafter  is  a  linear  function  of  the  applied  load; 

e 


aj ,  a9  =  constants 
P£  =  cracking  load 


aj  +  aj  P 


a.  +  a,  P  =  0 
1  2  c 


e  =  web  steel  strain 
v 


%  ■  a2  C-'c) 


(B.  21) 

V  L  V 

P  =  applied  load 

If  it  is  also  assumed  that  the  shear  force  at  a  section  is  proportional  to 
the  applied  loads,  then: 


C  =  constant 

V  =  shear  at  section 

V  *  cracking  shear  at  section 

f  =  web  steel  stress 
v 

E  =  linear  modulus  for  web  steel 
C'=  OE**constant 


e  =  C  (V-V  ) 
v  f  c' 

v 

ev  =  E 


f  =  C'  (V-V  ) 
v  cJ 


(B. 22) 


This  corresponds  to  the  same  form  accepted  by  the  ACI  code  (318-71)  and 
used  by  others  in  describing  test  results  [ B .  14 ,  17,  30,  38];  i.e.: 


b  =  width  of  member 


v  =  rf 


vy 


+  v. 


d  =  effective  depth 

A  =  total  web  ster'  area  at 
one  section  *  Afa  x  no. 
of  legs 


vy 


■  dfv(V'Vc) 


v 

c 


Vc 

bd 


(B.  23) 
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WEB  STEEL  STRAIN -€y 


*PC  *  LOAD  AT  CRACKING 


WEB  STEEL  STRAIN  x  I05-€V 


FIG.B.I6:  WEB  STEEL  STRAIN  AS  A  FUNCTION  OF  APPLIED  LOAD 
(DATA  FROM  REF.  B.  11,25]) 


1 


S  =  web  bar  spacing 

f  ,  V,  V  =  defined  above 
v  c 

If  it  is  assumed  that  the  crack  is  inclined  to  the  longitudinal  axis  at 

approximately  45°,  then  (d/S)  is  a  measure  of  the  number  of  bars  (n) 

crossing  the  crack;  i.e.,  n  =  d/S. 

According  to  observations,  all  of  the  web  bars  crossing  a  physical 

crack  are  affected  by  that  crack.  In  addition,  by  the  time  the  ultimate 

behavior  state  is  reached,  all  of  these  bars  can  yield,  [B.14,7]. 

All  of  the  bars  at  a  crack  should  be  included  in  the  total  effect, 

although  each  stress  cannot  be  predicted  independently.  Therefore,  the 

total  area,  (n.A^),  is  used  in  the  equation. 

The  average  stress  value  for  the  web  bars  crossing  a  crack  for  any 

value  of  shear  force  is  given  by  the  equation  for  fy  (equation  B.22) 

with  C'  =  S./dA  ;  i.e. : 

v 


f 

v 


S 

dA 


*  (V-V  ) 


v 

To  introduce  a  lower  bound  measurement  for  yielding,  the  data  plotted  in 
Fig.  B.17  is  used.  The  lower  limit  function  can  be  defined  by  the 
equation : 


vy 


-  *  .  (V-V  ) 

3  dA  1  V 


(B.24) 


which  is  valid  at  yielding  conditions  only.  To  introduce  f  at  any  stress 
state,  the  relationship  must  be  made  an  inequality;  i.e.: 


f  >  4  •  f  :  before  yielding 
vy  3  v  7  6 

4 

f  <_■=■*  f  :  at  or  after  yielding, 
vy  =  3  v  7  6 


125 


1 


FIG.  B  17  =  ULTIMATE  STRENGTH  MEASUREMENT  WITH  NORMAL 
WEB  REINFORCEMENT  (DATA  FROM  REF. [B.  30]  ). 
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(1)  Basic  f  lilure  criterion  :  based  on  an  equation  to  approximate  the 
web  steel  stress  at  a  crack: 


v  vy  v  J 

An  upper  limit  for  the  value  of  (V-Vc),  as  recommended  by  the  ACI  Code,  is 
not  included  in  the  criterion  for  two  reasons: 


(a)  Test  results  indicate  that  the  upper  limit  is  not  consistent  with 
physical  behavior  [B . 14 } ; 

(b)  A  natural  upper  limit  to  failure  for  large  values  of  (V-Vc)  is 
available  through  the  flexural  failure  criteria  . 

(2)  Axial  force  effect:  this  is  to  measure  the  effect  of  an  axial  force, 
tension  or  compression,  on  the  physical  behavior  of  web  reinforcement  after 
cracking.  Although  data  is  scarce,  reference  B.14  provides  results  to  indicate 
that  the  contribution  of  web  reinforcement  to  the  total  strength  is  approxi¬ 
mately  independent  of  both  the  axial  force  and  the  ratio  £a/u) .  These 
results  were  for  limited  va  i^tions  of  the  parameters,  but  they  form  the 
basis  for  the  criterion  statement.  The  measure  of  the  limit  state  for  the 
shear  strength  of  a  member  with  web  reinforcement  as  defined  above  is 
assumed  to  be  valid  for  a  moderate  axial  force,  either  compression  or 
tension.  It  should  be  mentioned  that  some  effect  is  already  built  into 
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the  expression  through  as  defined  by  V  in  section  B.3.2.1.  The 
quantity  that  is  assumed  to  remain  unaffected  is  (V-Vc)  or  (v*vcn)  •  The 
expression  which  includes  the  axial  force  effect  should  be  expressed  as 
follows : 


.  .  i 

3 

where  F  (x) 


Jr 

-r-*(F(x)  -  F  (x)J  <  0 
vy 


V  /  bd/f' 
cn  f 


(B.27) 


(3)  Loading  history  effects:  the  question  is  whether  unloading  or 
reversal  of  loading  has  any  effect  on  the  basic  behavior  of  web  reinforce¬ 
ment.  The  results  presented  in  reference  B.2  included  behavior  of  beams  with 
web  reinforcement  subjected  to  a  few  reversed  loading  cycles.  The  conclu¬ 
sion  stated  in  section  B.3.2.1,  part  3  is  applicable  to  this  part  also; 

i.e.,  the  repetition  of  a  few  reversed  loadings  does  not  cause  a  significant 
decrease  in  the  strength  of  beams  failing  in  shear  as  compared  to  monotonic 
loading  failure.  The  criterion  stated  remains  valid  with  load  reversal. 

(4)  Modification: 

(a)  Modification  by  coefficient:  the  web  strength  criterion  can  be 
altered  by  a  coefficient  which  changes  the  slope  of  the  function 
shown  in  Fig.  B.17;  i.e.,: 


f  =  C  •  -j—  •  (V-V  ) 
vy  5  dA  v  c 


The  suggested  relationship  developed  in  reference  B.14  to  fit  a 
specific  set  of  data  was  defined  to  be: 


1 . 75rf  =  v  -  v  ,  where  v  =  ultimate  shear  stress, 
vy  u  c  u 

This  equation  is  plotted  on  Fig.  B.17  and  it  appears  to  be  a 


suitable  upper  bound.  Accepting  this  function  as  defining  an 
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upper  bound  value  for  C$,  the  modified  form  of  the  criterion 
can  be  written  as: 


1. 


[F (x)  -  F*  (x)  ]  <  0 


for  (4/7  <  C5  <  4/3) 


(B . 28) 


and  is  included  in  its  modified  form  F*(x),  (equation  B.20). 

(b)  this  criterion  can  be  overridden  since  it  is  an  indirect  measure 
of  the  web  steel  stress.  If  it  is  ignored,  only  the  flexural 
behavior  limitations  are  in  effect. 

(5)  Assumptions: 

(a)  Web  bars  are  anchored  sufficiently  to  insure  yield  strength 
development ; 

(b)  Element  model  remains  valid  after  the  formation  of  tnc  diagonal 
crack  up  to  the  web  yield  state; 

(c)  The  web  steel  strain  is  linearly  related  to  the  applied  load  after 
crack  formation; 

(d)  All  bars  affected  by  the  crack  yield  at  the  measured  limit  state. 


B.3.3  AXIAL  FORCE  FAILURE 


The  limit  states  developed  in  this  failure  category  relate  to  the  dis¬ 
continuities  encountered  in  the  behavior  r"  a  member  dominated  by  normal 
stresses  caused  predominantly  by  an  axial  force.  The  distinction  between 
this  category  and  the  category  for  normal  stress  effects  due  to  flexure 
(B.3.1.)  is  that  the  point  of  zero  strain  at  a  specified  section  falls  out¬ 
side  of  the  physical  dimensions.  Hence,  the  state  of  stress  is  either 
tension  or  compression  over  the  entire  cross  section.  The  failure  modes 
which  require  a  failure  criteria  are  shown  in  Fig.  B.4. 
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The  failure  inodes  are  similar  to  those  for  flexure;  i.e.,  compressive 
failure  means  either  crushing  of  the  concrete  or  the  simultaneous  crushing 
and  bar  buckling  in  the  critical  compression  zone;  and  tensile  failure  is 
the  fracture  of  longitudinal  bars.  Therefore,  the  same  failure  criteria 
are  applicable.  In  the  axial  force  category,  the  secondary  stresses  are 
caused  by  bending.  But  since  both  are  included  in  the  normal  strain  state, 
these  secondary  effects  are  automatically  accounted  for  in  the  failure 
criteria.  This  is  similar  to  the  secondary  axial  force  effects  automati¬ 
cally  included  in  the  flexural  failure  criteria. 

B.3.3.1  CONCRETE  CRUSHING 

failure  criterion:  same  as  section  B.3.1.1. 

B.3.3.2  STEEL  FRACTURE 

failure  criterion:  same  as  section  B.3.1.2. 

B.3.3.3  BAR  BUCKLING  AND  CONCRETE  CRUSHING  SIMULTANEOUSLY. 

Failure  criterion:  same  as  section  B.3.1.3. 

B.3.3.4  STEEL  FRACTURE 

Failure  criterion:  same  as  section  B.3.1.  2. 

B.3.4  SUMMARY  OF  FAILURE  CRITERIA 

Basic  Form  of  Failure  Criteria  Expressions: 

[(Specific  Criterion  Value)  -(Computed  Model  Value)]  <  0 
All  criteria  are  expressed  in  a  dimensionless  form  except  for  B.3.1.1  as 
noted. 

( B.3.1)  Flexural  Failure 

(B.3.1.1)  Concrete  Crushing:  -  e£  <  0 

ec  -  Compression 
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'fl 


3  +  0.002f ' 

(r-~~ooo C)  =  °*0035 

c 


(1.  =  Cx  =  1.23) 
Default:  ■  1 

override:  not  possible 


Note:  Constants  in  the  expression  for 
* 

are  not  dimensionless,  even  though 
the  total  expression  is  dimensionless. 


(B.3.1.2) 

Steel  Fracture: 

ef 2  *  £s 

o 

VII 

e  -  tension 
s 

e^2  *  maximum  strain  value  defined  for 
stress-strain  function. 

override: 

not  possible 

(B.3.1.3) 

Concrete 

Bar  Buckling  and 

Crushing  Simultaneously: 

2 

r  .  ir  . 
2  16 

(f)2  .(^)  -  i.  i  o 

s 

fg  -  compression 

(l.  =  c2 

• 

VII 

Default: 

C2  *  2 

override: 

possible 

(B.3.1.4) 

Steel  Fracture: 

(Same  as 

B.3.1.2) 
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( tt«3.2)  Shea r-Flexure  Failure 


(B.  3 


(B.3.2.1)  Detection  of  the  Principal 
Diagonal  Crack  Leading  to 
Failure: 

F*(x)  -  F(x)  t  0 

F*(x)  =  {[1.5  +  C3  •  F2(x)  ]  *  ^7-.^2C3l}-  0  .  F^(x) 

(3.5  5  35. /6.) 

Default:  C3  *  3.5 


{Tension:  a  =  0.025(4.  -  C^) 
Compression:  ■  0.050(2.  +  C^) 


(0.  i  c4  5  1.) 

Default:  C^  ■  0 
override:  possible 


(B.3.2.2)  Detection  of  the  Principal 

Diagonal  Crack  Plus  Yielding 
of  the  Web  Reinforcement 
Leading  to  Failure: 


1. 


[F(x)  -  F*(x)]  i  0 


(4/7  *  C5  =  4/3) 


Default:  =  4 . / 3 


override:  possible 


3)  Axial  Force  Failure 


(B.3.3.1)  Concrete  Crushing:  (same  as  B.3.1.1) 
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(B.3.3.2)  Steel  Fracture: 

(B.3.3.3)  Bar  Buckling  and 
Concrete  Crushing 
Simultaneously: 

(B.3.3.4)  Steel  Fracture: 


(same  aa  B.3.1.2) 


(same  as  B.3.1.3) 
(same  as  B.3.1.2) 
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B.4  IMPLEMENTATION 


The  application  of  the  failure  criteria  is  made  at  an  equilibrium 
configuration  for  the  system  and  at  critical  sections  within  each  element. 
The  implementation  defines  the  critical  sections  used  and  describes  the 
general  application  procedure. 

B.4. 1  ASSUMPTIONS 

In  addition  to  the  assumptions  made  for  the  various  failure  criteria 
developed  in  B.3,  the  following  assumptions  are  made  with  respect  to  the 
implementation  of  these  criteria: 

1.  A  single  element  has  uniform  reinforcement  properties,  in  longi¬ 
tudinal  and  web  steel,  over  the  entire  length.  If  a  member  has 
variable  reinforcement  properties,  e.g.,  a  region  with  web  rein¬ 
forcement  and  a  region  without,  the  failure  criteria  will  be 
utilized  more  effectively  by  modeling  the  member  with  more  than 
one  element  to  reflect,  at  least  approximately,  the  actual  rein¬ 
forcement  uniformity. 

2.  General  beam  shear-flexure  behavior  can  be  effectively  related 
to  criteria  based  on  simple  beam  test  results  by  defining 
analogous  segments  within  the  model.  The  analogous  segments,  or 
shear  spans,  are  defined  by  lengths  measured  from  a  point  of  zero 
moment  to  a  nodal  point,  or  from  one  nodal  point  to  another. 
(Refer  to  Fig.  B.18). 

3.  The  most  probable  section  for  failure  in  an  element  is  defined  by 
the  critical  section.  A  critical  section  is  prescribed  for  each 
failure  criteria,  and  only  the  critical  sections  in  each  element 
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P|fP2«  CONCENTRATED  OR  EQUIVALENT  NOOAL  FORCES 
Q'b.C'd  »  NODAL  POINTS  FOR  THE  MEMBER 


FIG.  B.  18 :  POSSIBLE  FRAME  MEMBER  STATE  WITH 
SHEAR  SPANS  OEFINED. 
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are  checked  for  failure. 


4.  An  element  defined  with  a  length  leas  than  the  effective  depth 

d  will  not  be  checked  for  shear  failure,  but  will  be  checked  for 
flexural  and  axial  failure. 

5.  The  equivalent  distance  from  zero  moment  to  the  critical  section 
at  the  crack  location,  measured  by  (V(x)/M(x))  in  the  shear 
cracking  equation,  is  the  same  regardless  of  the  superimposed 
uniform  strain  caused  principally  by  an  axial  force. 

B.4.2  CRITICAL  SECTIONS 

Critical  sections  are  the  prescribed  section  locations  where  the  various 
criteria  are  applied.  It  represents,  in  each  case,  the  most  probable  lo¬ 
cation  for  failure  consistent  with  the  model  representation  of  the  actual 
member  and  applied  loading.  In  applying  the  criteria,  each  element  tested 
is  treated  as  an  independent  unit,  even  though  several  elements  may  be 
joined  together  to  represent  a  single  physical  member. 

Failure  due  to  flexural  and  axial  force  effects  are  related  to  maximum 
normal  strain  states  in  an  element.  For  the  element  model  used,  this  state 
occurs  at  one  of  the  end  sections.  The  end  sections  are  then  the  critical 
sections  for  these  two  criteria,  defined  in  B.3.1  and  B.3.3.  There  is  one 
exception  to  this  definition:  when  an  element  has  one  or  both  nodal  points 
located  at  a  junction  of  two  members  normal  to  each  other,  e.g.,a  horizontal 
member  intersecting  a  column  member,  the  critical  section  is  defined  at  the 
face  of  the  normal  member  rather  than  at  the  nodal  point.  This  can  be  seen 
in  Fig.  B.18  at  nodal  point  a  and  d. 

Failure  due  to  shear  effects,  with  criteria  defined  in  B.3.2,  is 
detected  by  an  indirect  measure  that  includes  material  properties, 


dimensions  of  Che  cross  section,  and  a  measure  of  shear  span  length.  The 
critical  section  Is  defined  by  a  conservative  estimate  of  the  location  of 
the  mid  depth  of  the  complete  diagonal  crack  within  the  shear  span  .  This 
location  is  specified  for  all  shear  span  lengths  in  Fig.  B.19. 

Others  have  suggested  locations  for  the  critical  section  in  a  shear 
span.  The  following  list  shows  a  reference  with  the  critical  location  used 
which  agreed  with  the  data  from  the  corresponding  tests: 


Reference 


B.10 
B.17 
B.27 
B.  30 


Critical  Section  Location 
X  (a  =  2d) 

_ cs  _ ' _ 

0.5  a 

0.6  a  (a  -  2d  for  a  >  5d) 
0.5  a 
a  -  d 


All  tests  involved  were  with  simple  beam  conditions,  and  the  cracking  stress 
equations  used  in  each  case  were  of  a  similar  basic  form  to  the  one  used  in 
the  criteria. 

Within  a  shear  span  associated  with  a  simple  beam  structure  with  a 
single  concentrated  force,  the  lowest  critical  cracking  stress  defined  by 
the  criteria  is  at  the  point  of  the  load.  That  section  is  the  most  con¬ 
servative  location  for  the  critical  section.  However,  it  is  not  realistic 
to  use  this  section  because  the  crack  has  to  have  some  room  to  form  before 
the  load  is  encountered.  Therefore,  the  critical  section  is  usually  chosen 
some  distance  from  the  point  of  the  load.  If  the  crack  is  assumed  to  have 
a  projected  length  equal  to  d,  then  (a  -  d)  is  the  most  conservative  realis¬ 
tic  choice.  Test  observations  indicate  that  the  location  is  actually  closer 
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SPAN-(l/d) 


(b) 


FIG.  B.  19 !  CRITICAL  SECTION  FOR  DIAGONAL  CRACKING 
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to  the  middle  o£  the  span  a.  The  location  uaed  with  a  specific  cracking 
equation  must  be  balanced  with  the  form  and  conservativeneas  of  the  equation 
itself  to  produce  realistic  and  conservative  results.  This  is  the  final 
decision  basis.  The  location  chosen  for  the  criteria  is  a  reasonable  balance 
between  the  most  conservative  location  and  the  observed  locations.  It  has 
also  been  checked  with  a  variety  of  test  data  and  compared  with  the  ACI 
Code  equation.  These  results  are  shewn  in  Fig.  B.20. 

A  summary  of  the  locations  for  critical  sections  associated  with  each 
failure  criteria  is  shown  below: 


Failure  Criteria 


1.  Flexural  failure: 


2.  Shear-flexure  failure: 


3.  Axial  force  failure: 


Critical  Sections  for  an  element 


End  sections,  at  the  nodal  points, 
or  at  the  face  of  an  intersecting 
member. 

Section  at  the  mid-depth  of  the 
diagonal  crack  -  defined  in  Fig. 
B.19. 


(same  as  for  flexural  failure) 


B.4.3  CHECKING  PROCEDURE 

The  two  general  element  states,  characterized  by  the  internal  stress 
distribution,  are  represented  by  element  (a  -  b)  and  element  (b  -  c)  in 
Fig.  B.18.  In  the  first,  there  is  a  zero  moment  condition  within  the 
element;  the  second  has  no  such  condition. 

In  the  first  element  (a  -  b),  the  flexural  and  axial  force  checks  can 
be  made  at  the  critical  sections  for  the  element  as  defined  previously. 
However,  the  shear  check  must  be  made  for  both  shear  spans  defined;  i.e. 
from  the  zero  moment  section  to  each  nodal  point  (shear  spans  1  and  2  in 


b. 
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(TEST) /»,;„( 


FIG.  B.  20 :  CHECK  OF  THE  DERIVED  EQUATION  FOR  SHEAR- 
FLEXURE  FAILURE  (DATA  FROM  REF.  [B.  14] ). 
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Che  figure).  Therefore,  the  zero  moment  eection  muet  be  located.  It  le 
defined  as  the  eection  with  a  uniform  strain  distribution. 

In  the  second  element  state,  (b  -  c) ,  the  flexural  and  axial  criteria 
are  applied  ae  before,  while  the  shear  criteria  can  be  applied  to  the  given 
state  directly.  The  application  of  the  shear  criteria  to  this  element  is 
equivalent  to  applying  these  expressions  to  a  shear  span  defined  by  the 
length  (g  -  c)  in  Fig.  B.18.  The  length  (g  -  c)  is  referred  to  as  the 
equivalent  shear  span  for  element  (b  -  c) ,  and  it  has  the  approximate  value 
given  by: 

Mb  ■  bending  moment  at  b 

M. 

V.  -  shear  in  (b  -  c)  L  -  - —  +  L. 

be  gc  Vbc  T>c 

L  “  length 

When  a  single  member  is  composed  of  several  elements,  the  question 
arises  concerning  the  adequacy  of  an  independent  check  of  each  element  to 
represent  the  behavior  of  the  whole  member.  In  other  words,  is  it  possible 
that  the  total  member  might  fall  if  it  Is  modeled  as  a  single  element, 
whereas  the  independent  checking  of  component  elements  might  not  detect 
failure?  It  has  been  demonstrated  by  physical  testing  of  restrained  beams 
that  the  most  likely  region  for  cracking  is  in  the  span  with  the  largest 
shear  force,  such  as  element  (c  -  d);  (reference  B.26  ). 

To  verify  that  the  criteria  expressions  would  be  consistent  with  this 
behavior,  an  idealized  analysis  was  performed  using  a  member  composed  of  two 
elements  subjected  to  the  possible  combinations  of  bending  moment  in  each 
span,  e.g.  the  two  spans  (a  -  b)  and  (b  -  c).  The  independent  checks 
produced  realistic  failure  predictions.  The  results  showed  that  the  cracking 
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is  less  likely  In  (b  -  c)  and  more  likely  in  an  adjacent  span  with  larger 

shear  force  (a  -  b)  as  the  momenta  at  b  and  c  approach  equality.  It  was 

also  shown  that  if  the  bending  moments  M^,  and  Mc  define  a  linear  function 

with  x,  and  M  *0,  then  the  failure  of  element  (b  -  c)  would  be  identical 
a 

to  the  failure  of  the  total  member  (a  -  c)  provided  the  same  critical  section 
is  used,  and  provided  both  spans  have  identical  properties.  Therefore, 
checking  each  element  of  a  member  Independently  will  accurately  reflect 
the  physical  behavior  of  the  total  member;  i.e.»  the  most  likely  failure 
region  will  be  detected  first. 

The  general  procedure  used  in  checking  for  failure  is  indicated  In  the 
flow  sequence  shown  In  Fig.  B.21,  where, 

critical  strain  values  -  strain  values  at  the  extremities  of  the 

end  sections  of  an  element; 

section  of  zero  moment  *  section  with  a  uniform  strain  state; 

equivalent  shear  span  -  distance  from  the  zero  moment  section 

(outside  of  the  element  length)  to  the 
section  of  maximum  moment. 

Flexural  failure  checks  (B.3.1)  are  required  at  both  ends  rather  than 
just  at  the  end  with  the  maximum  moment  since  an  element  may  be  reinforced 
unsymmetrically  with  respect  to  number,  diameter,  and  location  of  longi¬ 
tudinal  bars.  By  the  same  reasoning,  the  axial  force  failure  checks 
(B.3.3)  are  also  made  at  both  ends  of  an  element. 
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-B.21:  General  Sequence  for  Implementation  of  the  Failure  Criteria 
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NOTATION 


gross  cross  section  width 

effective  depth  of  longitudinal  reinforcing  bars 
total  depth  of  cross  section 

spacing  of  web  reinforcing  bars 

length  of  shear  span  for  checking  diagonal  cracking  failure 
shear  span  measured  from  zero  moment  to  a  concentrated  load 
reinforcing  bar  diameter 
radius  of  gyration  of  reinforcing  bar 

distance  to  diagonal  crack  intersection  with  tension  steel 
distance  to  critical  section  from  minimum  moment 

area  of  reinforcing  bars  as  a  group 
area  of  a  single  reinforcing  bar 

A^x  number  of  legs  ■  effective  area  of  web  reinforcement  bars 
(Closed  stirrups:  Av  *  2  x  A^) 

A  /bd  ■  Longitudinal  steel  pircentage 
8 

A  / bs  -  dimensionless  measure  Oi  web  reinforcement  area 
v 

d/s  ■  number  of  bars  crossing  a 'crack  (at  45°) 

unconfined  concrete  strain  at  0.20  f^  beyond  the  ultimate  strain 
unconfined  concrete  strain  at  0.50  f’  beyond  the  ultimate  strain 
strain  in  concrete  at  crushing  in  compression 
strain  in  concrete  at  crushing  with  modification  coefficient 


e  ■  strain  in  concrete  in  compression 

c 

ef2  ■  strain  in  longitudinal  steal  at  fracture  in  tension 

e  -  strain  in  longitudinal  steel  in  tension 

8 

ev  -  strain  in  web  reinforcing  bar 

f^  ■  ultimate  cylinder  strength  of  concrete 

✓F  a  mea  ure  of  the  tensile  strength  of  concrete 
f8  -  stress  in  longitudinal  reinforcing  bar 

f  -  compressive  stress  in  longitudinal  reinforcing  bar  at  critical 

condition 

fv  -  stress  in  web  reinforcing  bar 

f  -  stress  in  web  reinforcing  bar  at  yield  in  tension 

v  -  V/bd  -  nominal  measure  of  shear  stress 

vc  ■  Vc/bd  -  nominal  measure  of  cracking  shear  stress  without  axial 
force 

v  -  V  /bd  -  nominal  measure  of  cracking  shear  stress  with  axial 
cn  cn 

force  (tension  or  compression) 

Et  -  tangent  modulus  for  steel  stress-strain  response 

V(x)  ■  V  -  shear  force  at  section  x 

V  (x)  ■  V  -  cracking  shear  force  at  section  x  without  axial  force 
c  c 

vcn(x)  "  vcn  “  crackln8  shear  force  at  section  x  with  axial  force 

N  «  axial  force  in  an  element  (tension  or  compression) 

P  ■  measure  of  applied  load 

Pc  -  measure  of  applied  load  at  diagonal  cracking 
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M(x) 

F(x) 

Fx(x) 

F2(x) 

F3(x) 

F$(x) 


M  »  bending  moment:  at  section  x 

V(x)/bd/F 

Vc(x)/bd/F 

1000.  ]jd_  j  V(x) 

/f'*  |  M(x)  | 

Vcn(x)/bd^I 

V  (x)/bd/f"*"  “  measure  with  modification  coefficient 

cn  c 

N/bd/f7"  ■  nondimens ional  axial  force 
c 


Cr  C^,  C^,  C4,  »  modification  coefficients  (range  of  values  defined 

with  criteria) 


a  *  coefficient  for  axial  force  effect  on  shear-flexure  failure 
*  0.025(4.  -  C^) 

ac  -  0.050(2.  +  C4) 
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APPENDIX  C 
BIBLIOGRAPHY 

This  bibliography  contains  the  results  of  an  extensive  literature 
search.  The  purpose  Is  to  document  the  important  sources  of  information 
related  to  the  development  of  a  mathematical  model  to  study  the  complete 
behavior  of  a  reinforced  concrete  skeletal  structure  top  to  the  state  of 
collapse.  The  literature  associated  with  the  total  scope  of  work  is  so 
extensive  that  only  those  sources  related  to  the  development  of  the 
basic  model  are  included. 

The  organization  of  the  entries  is  by  basic  subject  division,  as 
shown  below. 

C.l  Material  Behavior 
C.1.1  Concrete 
C.l. 2  Bond  and  Anchorage 
C.l. 3  Reinforcing  Steel 

C.2  Element  Behavior 
C.2.1  Model 

C.2.2  Strength  Properties 
C.3  System  Behavior 
C.4  Solution  Process 
C.S  'Selected  Books 
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